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ABSTRACT 


This article is an account of the present position of theories of nuclear 
matter initiated some years ago mainly by Jastrow and Brueckner. Section | 
introduces the basic ideas. Section 2 discusses the variational approach. 
Section 3 develops methods involving partial summations in perturbation 
theory. The difficulties confronting the latter approach, and their partial 
resolution by ideas from superconductivity theory, are analysed in § 4. 

Much of the formalism is relevant also in the theory of individual nuclei 
(see for example the review by Eden 1959) and in fields other than nuclear 
physics (see for example The Many Body Problem, Les Houches, 1958). 
However, in the selection of material our criterion has been relevance for 
the ground state of nuclear matter. 
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$1. NuctEar MarTrER 
1.1. Properties of Nuclear Matter 


Ir has long been known that nuclei, unlike atoms, exhibit fairly well- 
defined radii and moreover that the nuclear volume is proportional to the 
number of nucleons. This led to the idea that, for some purposes, nuclei 
could be regarded as droplets of substantially incompressible ‘nuclear 
matter’. Ifthe coulomb forces, which prevent the occurrence of arbitrarily 
large nuclei in nature, are switched off, it becomes possible to contemplate 
in theory the existence of indefinitely large portions of this matter. The 
resulting spatial homogeneity greatly simplifies the theory. The pro- 
perties which such a theory should try to account for must of course be 
obtained by extrapolation from finite nuclei. 

The normal density of nuclear matter is obtained most directly from the 
electron scattering experiments (reviewed by Elton 1961 and Ravenhall 
1958). These are consistent with the assumption that all nuclei have the 
same density in a central unform region of varying extent, and, surrounding 
this, have a surface region of constant thickness (~2f). When so analysed 
(Elton 1958) the experiments yield for the central regions a volume per 
nucleon $779° with 

foe LIA fe ou in aed ca ee 


We refer sometimes to the ‘mean spacing’ of nucleons, 1-67), which is 
defined as the cube root of the volume per nucleon. 

The energy density of nuclear matter can be inferred from the Weizsacker 
semi-empirical mass formula (reviewed by Green 1958). This expresses 
the energy per nucleon of a finite system as the value appropriate in nuclear 
matter, plus corrections for finite size (surface energy), electrical forces 
(coulomb energy) and departure of the proton—neutron ratio from the 
value unity which we attribute to standard nuclear matter (symmetry 
energy). Simple estimates of each effect result in the formula 


(E/A) = —a, +a, A? +a; Z2A48 4.a,(N—Z)2A- . . (1.2) 


for the energy # of a nucleus with A nucleons, V neutrons and Z protons. 
Fitting the empirical data gives, according to Green, 


a,= 15-8 Mev, 
a, = 18-0 Mev, 
a,= 0:72 Mev, 
Q,= 24 Mev. 


(1.3) 


To account for the quantity a, is clearly one of the main aims of a theory of 
nuclear matter. To predict a, would require a discussion of finite nuclei, or 
at least of the semi-infinite system. The quantity ag in the coulomb energy 
is readily interpreted, in quantitative agreement with the electron scattering 
data, in terms of nuclear size. Finally a, is again a proprety of the infinite 
system, to be obtained by varying therein the neutron—proton ratio. 
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However, the treatment of symmetry energy in (1.2) is probably over- 
simplified. Cameron (1957) has proposed that the final term be replaced 
by 


(FSA) (aoe — das A) re Oo 


which includes a surface as well as a volume effect. According to Green 
the data then require 
Asy X 30 Mev, i 


Ass~40 Mev. | MS) 


Only dsy is relevant for nuclear matter. 

The analysis of electron scattering referred to above, and the Weizsacker 
formula, involve the idea of nuclear incompressibility. This of course is an 
idealization. A measure of the compressibility is given by the parameter 

ate ui E/A 
— Ser, y . . . . . e 1.6 

Yo dr? ( / ) 7 ( ) 
(where H/A is the energy per particle in nuclear matter and r, defines the 
density) evaluated at the value of r, for which 

@ (RJA) =0 (1.7) 

dr, ut . . . . . ° . . . 
(this just implies zero pressure, i.e. the absence of stress). Now the central 
regions of nuclei are not in fact quite free from stress; on the simplest 
picture the matter there is compressed by surface tension and distended by 
coulomb forces. Ona naive theory this results in a change of mean density 
given by 
202 
uf g Shoal Aga ioc Dae Re 
is Cr ea oO ee Kate 

away from the normal value defined by (1.7). For actual nuclei the second 
term is always the larger and 67, is negative. With K=170mev (the 
theoretical figure of Brueckner and Gammel 1958) (679/79) varies from 
about 0:02 in 1°Ag to about 0-01 in ’Au. This change would not be 
observable, but much smaller values of K would be ruled out by the 
electron scattering data. The rather smaller value 


Ki 70 MOVG A cae ht ot ean en ete? Le) 


was indeed obtained at one time from the ‘isotope shift’ in atomic spectra 
(Ford and Hill 1955) which can be attributed in part to the variation in 
mean density predicted by (1.8) for nuclei of equal Z and different A. 
However, Bodmer (1958) has observed that the unknown variation of 
symmetry energy with density is also important, and prevents the un- 
ambiguous derivation of a value for K. Other refinements of the theory 
are also important here (Breit 1958). Further estimates of compressibility 


214 J. S. Bell and E. J. Squires on the 


are obtained from semi-empirical theories of the nuclear surface energy and 
shape (Wilets 1958, Skyrme 1956) where K enters as an important para- 
meter. Comparing the theory with various versions of the Weizsacker 
formula, Wilets obtains values 


Kwx175, 218, 302 Mev. stem ve set ieee SANE 


The compressibility is therefore very poorly known. 

Once finite compressibility is admitted the possibility must be considered 
that the central regions of heavy nuclei are not in fact of uniform density, 
for coulomb forces will tend to produce a central depression. This is not 
decided empirically, for electron scattering is quite insensitive to central 
density (Ravenhall 1958, Elton 1961). In a very naive picture which 
ignores possible neutron—proton separation and forces due to density 
gradients, one obtains a coulomb effect at the centre of 2-5 times the mean 
value given by the first term in (1.8). With AK~170Mev this implies a 
central depression of density by some 20% in heavy nuclei. However, in 
the more elaborate theory of Wilets (1958) the effect is almost entirely 
suppressed, for reasons which are obscure. 


1.2. The Fermi Gas 


The simplest model of nuclear matter, and the point of departure of 
nearly all current theoriest, is the ideal fermi gas. This gives already 
first crude estimates of some of the quantities mentioned above. Consider 
a cubical box of volume Q, which is finally made very large, on the boun- 
daries of which wave functions satisfy the conventional periodic boundary 
conditions. Itis assumed that for sufficiently large Q the precise boundary 
conditions are unimportant{, and the periodicity condition is convenient 
because it leaves the theory accurately translationally invariant. The 
permitted states of individual particles are plane waves with momentum 
vectors 


(k,, ky, k,) =2n(1, m, n)\Q-¥8 


where J, m, n are positive or negative integers. The ground state is a 
determinantal wave function in which all plane waves with |k | less than or 
equal to the fermi momentum &, are occupied, each with four particles— 
corresponding to the four states of spin (including isobaric spin). It is 
useful to denote occupied states by symbols m,m’, .. , unoccupied states by 
I,!’..., and states in general by k, k’.... Moreover it is usually con- 
venient to suppose that these symbols specify spins as well as momenta, 
to avoid encumbering formulae with special spin indices. 


RARE EL eee SS RE OLN OMe eA ee 

+ See however Trainor et al. (1958) and Trainor and Sood (1959). 

{It would be desirable to start with a theory of large finite nuclei and 
demonstrate that the periodic boundary condition model correctly represents 
their central regions. This has not been done. Some preliminary observations 
are made by Eden (1958). - 
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The fermi momentum is related to the density parameter 7, through the 
mean volume per nucleon 


ST gi oP bake 
For large Q there are if 
Qdm (27)-3 Gartmciatn cmt a ese UT) 


states available in a region dm of momentum space, where 
dm=dm,dm,,dm,. 


Then, with a factor 4 from spin summation, 


or 
kept = (977/8)'3 = 1-52 (1.12) 
With ro=1-12f, 
Hope fal celica cont ct ogg tie goss (Ls 08) 
The corresponding kinetic energy is 
je Po Mie 38 Mev le 8) ek ae vk v.47 41,14) 


and the mean kinetic energy per particle is - 
BU 2M eS MOV ort oy, Boys) 322 (V518) 


In the absence of interactions the energy per particle does not of course 
exhibit a minimum at any value of 7). But with r, fixed at the desired’ 
value, values can be obtained for compressibility and symmetry energy. 
The former is 

K =1,2(d?/dr,.?)(3k 7/101) = (18k,"/10M) 


Col Sa MOV eile bas en ee aad wR Cra ea yrees a) oC dO) 


The symmetry energy, obtained by considering slightly different fermi 
momenta for neutrons and protons subject to fixed total density, is 


Osy=(ke2/6M)wi3Mev. . ... . . « (117) 


These are already of the right order of magnitude, but comparable contri- 
butions arise from internucleon forces. 


1.3. Interactions 


The most fundamental procedure at this stage would be to let the fermi 
gas interact with a meson field and to evaluate the resultant energy shift. 
Such a programme is hardly practical, in view of the very limited success 
of meson theory even in the two-body problem. Nevertheless some 
speculative attempts have been made, based on the extreme simplification 
that the discrete nature of the nucleons as sources of the meson field is 
ignored and that they be replaced by an averaged density distribution 
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(Skyrme 1955, Johnson and Teller 1955, Schiff 1952, 1954). There results 
a collective meson field, and individual nucleons interact with this rather 
than with one another. This picture of bulk nuclear matter is not 
incompatible with the action of residual two-body forces between particles 
outside closed shells, as required in the shell model when it goes beyond the 
extreme single particle version. Nor is it incompatible with strong two- 
body correlations in the lower density surface region of a nucleus, which 
might possibly be enough to account for various high energy pheno- 
mena which have been put forward as establishing the existence of such 
correlations in nuclear matter (Brueckner et al. 1955). However, while 
such theories are not ruled out by experiment, and certainly not by meson 
theory in its present state, they are perhaps rather implausible. The 
mean spacing of particles in nuclear matter is slightly greater than the 
pion compton wavelength ; thus, while the meson clouds of several nucleons 
will often overlap. it seems unlikely that they fuse into a collective whole. 
Whatever the case may be, these theories have developed little recently, and 
we do not consider them here. 

The more traditional approach has been to suppose the meson field 
formally eliminated, leaving an effective hamiltonian involving nucleons 
only, and to determine the latter as well as possible by experiment. There 
will be terms involving the coordinates of two, three, four, ... and arbitrarily 
many nucleons. Much can be learned about the two-body interaction 
from nucleon-nucleon scattering experiments and deuteron properties. 
Something could then, in principle, be learned about three-body forces from 
3He and °H and from nucleon—deuteron scattering, but this would require 
rather well-determined two-body forces and a reliable solution of the 
wave-mechanical three-body problem (see for example Blatt and Derrick 
1960)+. In fact it seems impossible to say more about the experimental 
existence of three-body forces, and of many-body forces in general, than 
that so far no clear evidence for them has been recognized. Guidance 
from meson theory is equally indecisive. Drell and Huang (1953), where 
references to earlier literature are given, considered the problem in lowest 
order perturbation theory with the relativistic PS-PS meson theory. 
They found the three-body force to be of very great importance in nuclear 
matter. However their leading terms involved s-wave meson—nucleon 
interaction, for which perturbation theory is now known to be grossly 
wrong. In this respect the work of Fujita and Miyazawa (1957) is more 
realistic}. They use no-recoil meson theory and by dispersion relations 
relate the simplest contribution to the three-body force to observed meson— 
nucleon scattering. The expectation value of this force in the ideal fermi 
gas at normal density is found to be only 0:2Mey. As the authors say, the 
of OS ee Sa ie ris SESS Ratti REA I at 

} These authors now find that the Brueckner-Gammel potential of this 
section, unsupplemented by any three-body force, is in definite disagreement 
with empirical properties of the three-nucleon system (Derrick et al. 1961). 

{ See also Brueckner et al. (1954) and Smith and Sharp (1960). The latter 
authors abtain a contribution to 4/H of just over 1 Mev at normal density. 
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result is not to be taken too seriously. The theory is reliable for at most 
the long-range part of the force and probably, as for the two-body force, 
quite inadequate for the bulk of it. The calculation does however provide 
some excuse for the customary position of leaving till later the considera- 
tion of many-body forces in nuclear theory. This position we adopt. 

It remains to make some brief remarks about the two-body potential 
(reviewed by Phillips 1959, Gammel and Thaler 1960 and MacGregor et al. 
1960). This is not determined observationally, even in principle, but 
involves an element of convention (see Phillips 1959 and the following 
paragraph). It is usually taken to be velocity independent, or nearly so, 
and simple assumptions made about its shape. Such a potential can still 
contain numerous parameters; the present data require the introduction 
of many, without being able to fix their values very definitely. However, 
some features appear to be clear. At least in the singlet s-state the inter- 
action appears to be strongly repulsive at short distances. This is usually 
represented by making the potential infinity positive within a radius of 
about 0-4f—the so-called ‘hard core’. Outside this there is an attractive 
central potential mainly in even states (a Serber force). The quadrupole 
. moment of the deuteron requires a tensor force, and in order to fit both 
high and low energy scattering it is probably necessary to add an interaction 
between spin and orbital angular momentum (this is a departure from 
velocity independence). In nuclear matter the tensor and spin orbit forces 
do not contribute to the energy in first order, and are often omitted in 
preliminary studies. So far only Brueckner and Gammel (1958) have used 
a potential with all the features listed}. This they obtained by substantial 
modification of one due to Gammel and Thaler (1957) which did not fit low 
energy data (the effective range region). Their potential is a hard core 
(in all states) of radius 0-4f outside of which 


v= (4—4P° P’){ — 877-4 (4 + 4P”) exp (—2-091r)/2-091r 
—434 (4—4P°)exp(—1-45r)/1-45r 
— 159-4 Sexp (—1-045r)/1-045r 
— 5000 S’ exp (—3-7r)}/3-7r 

+ (4+ 4P? P*){ — 14-0 (+ 3P’) exp (—1)/r 

+130 (4—4P’) exp (—7r)/r 
+22 Sexp (—0-8r)/0°8r 
— 7315 S’ exp (—3:7r)}/3°7r 


in Mev, with r in fermis, where 


S=3 


O,.1G,.r 
14 — 6, . 6s, S!= 3(hy— 2) X (Pi — P2)- 


r 
The precision with which the parameters are fixed by experiment is not 


stated. However the weak triplet-odd interaction in particular (the 5th 
ee ne ae ale ae ee eee 


+ Also Scott and Moszkowski (see § 3.5). 
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term above) is fixed quite poorly. It is heavily weighted in the nuclear 
matter calculation, and the quoted value was chosen to obtain there the 
desired energy per particle, which would otherwise be uncertain by several 
mev (J. Gammel, private communication). The Brueckner-Gammel 
potential does not at long range agree with the one-pion exchange potential 
of meson theory. This asymptotic part of the mesonic potential is now 
thought to be both theoretically reliable and experimentally supported 
(MacGregor et al. 1960). It is to be hoped that detailed potentials with the 
appropriate long-range character will be forthcoming and the nuclear 
matter calculations repeated. 

The hard core, as will be seen, is a considerable complication in the many- 
body problem. It has occasionally been suggested therefore that some 
velocity dependent (or non-local) potential be used instead (see, for 
example, Peierls 1960), which could equally well become repulsive at high 
energy. The empirical success of the hard core does not exclude this - 
possibility, for from any potential which fits the data another is generated 
by any canonical transformation which leaves unaltered the asymptotic 
wave function. It should be noted however that such a transformation 
alters the interpretation of the wavefunction and in general changes the 
expressions for quantities like, for example, the deuteron quadrupole 
moment. Moreover, potentials equivalent in the two-body system are not 
immediately so in the many-body case. This is compensated by the fact 
that the transformation, generalized to the many-body system, alters also 
the many-body forces or introduces some if they were absent. The 
three-body force of Fujita and Miyazawa, for example, is complementary 
only to the static two-body force (inevitable in static meson theory) 
explicity assumed in their work. However, these matters have received 
little attention. The only non-static potential we will refer to is the so- 
called ‘separable’ potential, and that less for its ability to fit data 
(Yamaguchi 1954) than for its great simplicity in formal investigations 
like those of § 4 where realism is a secondary consideration. 


1.4. Perturbation Theory 


Suppose for the moment that the nucleon interaction is a non-singular 
two-body force of the type for long conventional in nuclear physics. A 
non-singular potential is here one which has a finite fourier transform ; this 
does not exclude for example the Yukawa shape. Such a potential can 
fit data in the region of effective range theory, and some features of higher 
energy scattering can be obtained with the Serber exchange mixture. 
Perturbation theory can be used to estimate the resultant energy shift in 
nuclear matter away from that for the ideal gas. This was carried to second 
order by Kuler (1937) for a gaussian potential, by Huby (1949) for a Yukawa, 
and by Thouless (1957) for an exponential well acting in S-states only ; 
their general results are similar. It might be feared that the series would 
converge slowly if at all, for nuclear forces are not weak. In nuclear matter 
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particles collide with relative velocities up to (2k,/M), corresponding to 
incident energies in scattering ge up to 


> (ae) =4 20 = 150mev. Rtg tamer eteas | © Bet Re) 


Born approximation in collision theory is poor over all of this range (Jost 
and Pais 1951). Despite this the first approximation to the energy shift in 
nuclear matter, as judged by the size of the next term, is reasonable 
(Swiatecki 1956, Bethe 1956). Thus the perturbation results, and it is their 
main interest, support the idea of Weisskopf (1950) that the ideal gas is a 
much better starting point than might have been expected because the 
exclusion principle suppresses much wave function admixture that could 
otherwise occur. 

In first order the energy shift is the expectation value of the potential in 
the unperturbed state; i.e. 

(AZ),=4 > {<{mm’ |v] mm’)— (mm’ |v|m’m)} (1.19) 
in terms of matrix elements of the interaction, v, between pairs of single 
particle states normalized in the box. The presence of the exchange term 
will often be indicated by the notation 


(AZ),=4 > <mm’|olmm’—m'm). . . . (1.20) 


In calculating Si sums the form 


AE),=4 2, {mm |v(1—P’P?P")|mm’) ce Rael) 


is useful, where (P’, P’, P’) are two-body exchange operators in space, spin 
and isobaric spin. Consider for simplicity a mixture of only Wigner and 
Majorana forces 
v=(W-—MP?’ P*)u, 
Wat, I 


where u depends only on interparticle distance. The spin sum implied 
‘in (21) gives a factor which is the trace over the spin states of two particles 


trace (W —MP?P’)(1—P°P"P’) = 16{(W—4M) + (M— ZW) P. (1.23) 
The matrix element (k,k, |v|k,’k,’) has the spin independent factor 


(1.22) 


0] dey drgu(r 1—"o) exp {i(k,’ — k,) . ry +7 (Ky — ky)rg} 
= 0-2 8(k, +k, —k,’ —k,’) @(k,—k,’), . . . (1.24) 
where 
i(k) = | dea (7) exp (ihe). toe eee Lpae) 


Converting momentum sums into integrals by (1.11), and with 


dm 
(2) 


A=40 (1.26) 
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we obtain for the energy shift per particle 


AE dm dm’ a ke ; dm 
(OE), =2 [aap eegp 7 0 0) + 417) atom — mt | ES. 
(1.27) 


At high density (large k,) the exchange term, involving @(m—m’), 
becomes relatively negligible because of the fall off of a(k) for large k. 
The direct term, involving (0), is simply proportional to k,?, and at high 
density dominates also the kinetic energy (ock,*). Since @(0) is negative for 
forces attractive in S-states, it is then possible by going to arbitrarily high 
density to obtain an indefinitely low energy per particle, and the system is 
unstable against collapse, unless 


(W=tM S05 Bas, Cer eas 


This condition is quite accurately necessary, because the result of first 
order perturbation theory (not of higher orders) has a variational signi- 
ficance; the expression for the total energy is just the expectation value 
of the total hamiltonian (in the unperturbed state) and the true energy is 
lower. It can be shown to be sufficient as well as necessary (Wigner 1936). 
For a general, account of such ‘saturation conditions’ see Blatt and 
Weisskopf (1952). The Serber force (W = M=34) does not satisfy (1.28), 
and nuclear saturation is now attributed only in part to the exchange 
component in the force, which is supplemented by the hard core. 

To obtain a numerical value for (AH/A), at arbitrarily prescribed density 
we will use a Yukawa potential 


exp(—p7r) ¢ 4nV 
u=V ——_—_.,, k)= ~—,.. .. . (1.29) 
pr Ee pe + pk? ae 
Suitable parameters (used for example by Bethe 1956 and Brueckner 
1958) are 
V=49mev, pt=1-19f. . . . . . (1.80) 
These fit low energy singlet-state scattering. Performing the integrals in 
(1.27) (Huby 1949) we get 


AE BP dees ; ; 
(2), 3 f-unn) or- ai (2) 
with (1.31) 


8a 
g* 


1 1 1 
Q,(a) = ( 2a— — pag bs So 2)__8tan-19 
1(@) ( sa) + (; + iam) log (1 + 4a?) — § tan—! 2a. 
For ro=1-12f and (ky/u)=1-61, Qox3-7 and Q,x 1-2. 
The Serber mixture then giver 

(AH/A),~ —48mev. 2. 2. 2. 4.00% (1.32) 
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With the kinetic energy 23 Mev per particle, the binding energy per particle 
is 20Mev. For the just saturating mixture W=1, M= #, we would have 
obtained —2mev. The experimental value is 16 Mey. | ; 

Consider now the second order energy shift. This results from a two 
stage process; first a pair of particles is excited into states which, because 
of the exclusion principle, may be taken to be above the Fermi surface; 
these particles must then be returned to their original states, possibly 
interchanged. Thus 

(AB eT Ss SNe a Se ao mim ) 4 
Lea ee (m? + m'2 —[?—1')/2M 

The sum is over all spin states; |/l’)<Il' | is just the unit operator in spin 
space, so the required spin summation is 


4s trace (W — MP P*)? (1— P? Pt P’) 
=(W?+ M?-4WM)+(2WM—4W2-4M?)Pr. 
=a+ BP". 


(1.33) 


Using momentum conservation, implied by the 5-function in (1.24), we 
introduce 


k=Il—m=m'_-I’ 
ome ee i (1.35) 


k’=l—-m’=m—-l’=k+m—-mM’. 
Then, converting sums to integrals, 


(a8) -- har | dm dm’ dk odi(k)ii(k) + Ba(k)a(k’) / | dm 


A (ar)? (2ar)? (277)8 2k .k’ 
im+kl >ke 


Im’—kl>ke aoe (LO) 


where, as always, m m’ integrations run only below the fermi surface. For 
the Yukawa force 


(F) = bo) san(t)} om 


The integral Q, has been reduced by Huby to a lengthy combination of one 
dimensional integrals involving elementary functions. We refer to his 
equation (26) for the expression; his diagram includes a plot of (0-52Q,) 
against 7p (79 in his notation). The exchange integral Q, is more difficult 
but can be bounded by observing that (AZ), must be negative for arbitrary 
exchange mixture (because the first order shift is an upper bound) and that 
di(k) (1.29) has everywhere the same sign ; thus 


0<.0,<203 . (1:88)e) 
For the above parameters and at normal density Q,+ 0-50. Then, with 


the Serber mixture, 
4-6mMev >(AH/A),>2:3Mev. . . . . (1.39) 


This is only between 5% and 10% of the first-order contribution. At high 
density ©, increases only logarithmically with k,, while Q, goes as k,*, so 
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the convergence improves. To some extent this is to be expected simply 
because the two-body collisions are increasingly energetic. However, the 
exclusion principle plays a large part. To see this one can drop the 
restriction on the / integration in (1.36), taking the principal part of the 
resulting singular integral; , is replaced by 

(GQ a)= 7 4(3- a) Ok (5 + =a) log (1 + 4a?) — tan“ 2a ; 


Fabia gk OE): 


This increases more rapidly (k,) for large ky. Moreover, at normal density, 
with a= 1-61, Q,’ is between four and five times at large as Q,. The second- 
order theory has been extended to tensor forces by Levinger et al. (1960), 
who consider several possible shapes. The tensor force does not contribute 
in first order, vanishing when the spin average is performed. In the same 
way there are no cross terms in second order between central and tensor 
forces. However, the pure tensor contribution can in second order be very 
substantial. Ifa ‘weak’ tensor potential is used, i.e. one just sufficient to 
account for the observed stronger interaction in triplet than singlet states 
in the effective range region (assuming equal central forces), it gives an 
energy shift of about — 5 Mev per particle. The low energy data, including 
the deuteron quadrupole moment, can equally be fitted by a ‘strong’ 
tensor force and a much weaker triplet central potential (with about 0-6 of 
the singlet strength) ; the tensor shift in second order is then about — 10 Mev 
per particle. By contrast if no tensor force is allowed but the triplet 
central force increased appropriately (by a factor of about 1-4) the second- 
order shift per nucleon increases by only about—1 Mev. The first-order shift 
of course increases also, by about — 9 Mev per particle. Thus tensor forces 
tend to dominate in the second-order term, and appreciably worsen the 
convergence. 

Some discussion of third-order perturbation theory is contained in 
§3.6. With central forces, according to de Dominicis and Martin (1957 a), 
the rate of convergence is maintained. The effect of tensor forces in third 
order terms has not been discussed. 


1.5. The Hard Core in the Many-body Problem 


Ordinary perturbation theory is inapplicable for realistic nuclear 
potentials because of the hard core; even if this were taken to be finite it 
would have to be so large that the convergence of the perturbation series 
would be very poor. Before introducing the two principal methods which 
have been used to handle this situation we note that, since the radius of the 
core (~0-4f) is small compared to the mean separation (1-6ry%1-8f) of 
nucleons in the unperturbed fermi gas state, the effect of the core may, in 
some sense, be regarded as ‘weak’. Thus we expect that, in spite of the 
strong potential in the core region, the system retains many of the features 
of the unperturbed state and, for example, remains more closely akin to a 
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gas than to a liquid. This will be illustrated in §2 where it is explicitly 
demonstrated that a hard core of very small radius has negligible influence 
on the properties of the system. There is some experimental support 
that this situation holds to a large extent in the nucleus from the fact that 
the mean free path of low energy nucleons incident on nuclei is known to be 
large compared to the mean separation. Thus the unperturbed fermi gas 
is a reasonable starting point. 

The first method is based on the use of trial wave functions and the 
variation principle. As has been stated, the results of first-order perturba- 
tion theory can equally well be regarded as arising from the use of a deter- 
minantal trial wave function of plane wave states. With a hard core 
potential this only gives an upper limit of (+ 00) on the value of the energy 
per particle. By relaxing the translational symmetry of the single particle 
states it is of course possible to obtain a more reasonable bound}—the 
individual nucleons would have to remain in spatially separated states 
so that their hard cores did not overlap. The considerations above would 
suggest that such a highly correlated wave function would be very unreal- 
istic. An alternative is to multiply the original determinantal wave func- 
tion by a factor which vanishes whenever any two particles approach within 
the hard core radius. The simplest form of this correlation factor is a 
product 


ae eee eee et ree Seni Cr (Tea) 


of factors for each interparticle distance, where f vanishes over the appro- 
priate range. Outside this range f can be determined in principle by 
minimizing the energy. ‘This is the method proposed by Jastrow (1955). 
It is not without difficulties in practice. Itis the main subject of § 2. 

The second method is a development of the perturbation rather than the 
variation approach. It can be introduced by noting that the second-order 
perturbation energy can be incorporated in the first order by replacing there 
the matrix elements of v by those of the operator 


v |< |v ie 
‘+ 2e-PpM—-P pM eed 

with €=(m?+m")/2M. This can be written more compactly 
v+v woe SIN SMI sane eI EeO) 


where 7’ is the kinetic energy operator and the two-body operator projects 


+ That the Hartree-Fock state with least energy need not have translational 
symmetry has been stressed recently by Overhauser (1960). Our example shows 
that the lowest energy Hartree-Fock state is not necessarily then the best 
starting point for the theory. See also Kohn and Nettel (1960) and Broutt 
(1960), who suggest that Overhauser’s specific proposals may not be relevant 
in three dimensions, and Henley and Willets (1961), Sawada and Fukuda (1961) 
and Bloch (19600 a). ; 
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onto states above the fermi surface. Ignoring the exclusion principle as 
in obtaining (1-40) we would have instead 


Osh al Ba a ee 


where P indicates the principal part. Now these are just the first terms 
in the expansion of the ‘ reaction matrix’ of collision theory, closely related 
to the scattering operator (for a convenient introduction to these ideas see 
Brueckner 1958). The complete series is generated by the integral equation 


lg 0 
PRUE)... +s (1-45) 


70 / oe 
KE) vhs 


The operator K® remains quite finite and well defined even when the 
potential has a hard core (the matrix elements are linear combinations of 
the tangents of the scattering phase shifts) although each term in its 
perturbation series becomes nonsense. Now when higher order terms in 
perturbation theory for nuclear matter are examined a sub-set is found 
which is generated by the equation 
@ 

K(@)=v+v ep ke). oNiw ot Soh he, ee MERE 
This is very like (1.45), and it can be hoped that there also the solution 
remains sensible for singular potentials. It represents roughly speaking 
the effect of treating two-body collisions in all orders of perturbation theory 
while ignoring many-body effects. The operator X, rather than », is then a 
reasonable measure of the effective strength of the interaction, and it is 
reasonable to rearrange the whole series in terms of it. Of course such 
rearrangements of series with unexamined convergence properties are not 
mathematically respectable. We take the usual physicist’s view that 
they are permissible until things go manifestly wrong—which they do. 
This approach, due mainly to Brueckner, is developed in §3. What goes 
wrong is discussed in § 4. 


§ 2. TRIAL WAVEFUNCTION MrrHops 
2.1. T'wo-body Correlations and the Cluster Expansion 


Here we take up in detail the method of correlation factors mentioned 

briefly at the end of the last section. We introduce the trial wave function 

Y= I f(r) ® 

i>j 

=FQ, Soe Se Tat Loree, sb 

where 7; =|r;—r; |, © is the fermi gas state, a determinant of plane waves as 

discussed in $1.2, and / is a function which is real (using time reversal 
invariance) and which satisfies 

Tit) =x Oe A pe ge 
1, (ft Sey ° . e . e ° oT (2.2) 
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The last condition is imposed to justify the neglect of surface integrals in 
partial integrations used to obtain the Euler equation. Subsequently we 
allow F to be arbitrarily large. The form for f given in (2.2) is appropriate 
for a nuclear potential containing a hard core with radius c in all states, as is 
usually assumed. We shall frequently denote f (7;;) by f (¥). 

We require to determine f(r) by minimizing A, given by 


___ [(F®)*H FO 
is | (FO)* FO 


with respect to variations in f(r). In (2.3), and in future equations, we 
use the notation that, when the integration variables are not indicated 
explicitly, an integral sign stands for integration over all nuclear coordinates 
and summation over all nuclear spins. H is the complete hamiltonian for 
the nucleus. Unfortunately it is not possible to evaluate Af in closed form 
and it is necessary to use some approximate method. This means that the 
variational principle, Z < H, is only approximately valid. It is important, 
therefore, that f should be varied within some restricted class of functions 
for which the approximate method of evaluating H is adequate, and, as we 
shall see, the omission of such a restriction can cause trouble. It is by no 
means certain that the restricted class is sufficiently large to allow a good 
estimate of H. 

The approximation used to evaluate Al is based on the expectation that 
the two-particle correlation, caused mainly by the hard core, disappears 
rapidly outside the range of the core, i.e. [f?(7)—1]| tends rapidly to zero 
for r>r,. This enables one to make an expansion in the range of the 
correlation function divided by the mean separation. Such an expansion, 
the ‘cluster expansion’, is familiar in the classical theory of imperfect 
gases. The method was first used in a quantum mechanical variational 
problem by Dingle (1949), who applied it to a hard sphere boson gas. 
Usually only the lowest order term is kept, which means that simultaneous 
correlation between three or more particles is ignored. Some indication 
that this might be a useful approximation is given by the fact that 
(c/7) = 0-2, so that, if f is taken to be zero inside the core and unity outside 
it, the higher terms of the cluster expansion should be small. Actually 
such an fis very unrealistic and we shall see that more reasonable forms give 
significant corrections. 

The equation for H (2.3) can be written 


—1iM Sor V2(F®) 3 & | oF (ij) FO 
u + v9 


florr flor 


where »(ij) is the interaction between nucleons i andj. To consider the 
potential energy term it is convenient to write v(t) in the form 
+ 


v(ij)= & Un (tg)Onl) - » 2 + + © (25) 


n= 


(2.3) 


A= (2.4) 


P.M.S- S 
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where 0,=1, O,=6;.6;, O;=7;.t; and O,=6,;.0;%;.7;. We do not 
include tensor forces since, with the wavenine Onn (2. 1), they will not 
contribute to the energy. Then we can write for the potential term in 
(2.4) 


pA | Yo,(r)G f(r) dr, 
where 
AD acl 2 ‘3 2772 ‘ 2. 
Gul) = = ay | PP On(12) 0 ee (2.6) 
an 


= (A/). 
We use the equations 


| o*F V2(FO) =} i] [FAO* V20+ 0 V20*}+ 2/0 PFVZF + Vj 2 VF2], 
Vi f(y)? 
VFi= YS pe. 
df ~ f(y? 


and 


i f(y) 


to write the kinetic energy term in (2.4) as 


vie SO es. yws.0(2), 


[47. - fo fare V2f)G,+4f Vf. VG, +4f Vf-H 
where 7')( = 3k,?/10M/) is the unperturbed kinetic energy per particle and 
“as Wie 0 : SNe 
He)= Saad) lOP/[F IOP 27) 


Integrating by as and adding the potential energy we have finally 
Aj/A=7,+ 2 siz | arte (Vf)?+4(VG,—H). Vi}+ 2S lars, Ge 


ott eecta eee 


Provided the rangeof fis sufficiently small it will be a good approximation 
to replace F? and F?/f? by unity in the equations for G and H. 


Then, H= VQ, 
and Gy =a (kr) Ve eee iS 4 eS «4 (2.9) 
where d(x) =1—412(x), 

d(x) = d3(x) =4d,4(~) = —312(x), } e e e (2.10) 
with U(x) = 32-8 (sinw—mcosz) . .« 0 «© 6 (2.11) 


The first approximation to (2.8) is then 


AlA=Tyt pdr) (WiP+4 Yom daf*t, 5 (2.12) 
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This is the first term in the general cluster expansion of Iwamoto and 
Yamada (1957a). A different form is given by Hartogh and Tolhoek 
(1958) and by Aviles (1958), who develop a systematic expansion in powers 
of the density parameter k,. To obtain the first term in this the functions d 
in (2.10) are replaced bys their values at the origin kyr=0. The latter 
expansion converges very poorly for nuclear matter (Temkin ue because 
the force range is comparable with the internucleon spacing ky 

As a trivial application consider the pure hard core Hotsatial Take a 
correlation function f (r/c) which is zero for argument less than unity. The 
potential energy term vanishes exactly, and (2.12) gives 


(AJA)=7,+ Pauld) or | dex? { f'(e 


If the cluster expansion converges this is an upper limit for small c. For this 
purely repulsive potential 7’) is a lower limit. So we confirm formally that 
a sufficiently small hard core has negligible effect on the energy. 

In the general case Edwards (1960, private communication) has found 
that the cluster development can be avoided and (H/A) evaluated exactly 
if the correlation factor IIf? is somewhat modified. This has the advantage 
of preserving the variational significance of the method. Unfortunately 
not enough correlation is retained to give finite results with hard core 
potentials. 


2.2. The Euler Equations for the Correlation Function 


Consider now the minimization of H/A, given by (2.12), for arbitrary 
variation of f. Taking for simplicity a spin independent force, we must 


“minimize | 
VI pg yf 2.13 
Ty+pldrdyJ—o +hofth. . 2... (2.18) 


The Euler equation is 


va (deN et iaa 0,0. s ora et ah ons 7 (2.04) 
which is to be solved with f(R)=1. When the restricting radius & is made 
arbitrarily large the solution has pelir ik form 

f(r)-1~ -a/r Se i tauie et Monee) 


where « is the ‘scattering length’ for the modified Schrodinger eqn. (2.14). 
Thus f has a long-range part which vitiates the cluster expansion. Never- 
theless, (2.12) can still be justified as a consistent first approximation if « 
is small. For then a weak non-singular potential v(r) can be chosen in 
many ways such that 
v=v' +6v , 

where v’ has zero scattering length. It is not hard to see that (2.12) is 
correct to lowest order in the range of v’ and in the strength of dv. 


$2 
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The energy shift per particle can be given in terms of a. By partial 
integration in (2.13), using (2.14), 


BlA=p[arv {obs Ve} 


oe ea yet! 
= psy Let ( (d, frp’) 


=p (2rafM).: inc as Se Se eer Lay 


For the pure hard core (2.14) gives 


pace Pai Oe 3 
t = | ret ik Hi (ba) ee 


For large 7, to first order inc, 


f=1—d,(0) (¢/r) 

whence 

2p 3c 
ng PEs RS na ae ee 
peels (2.18) 
plus terms of higher order in c. For small ¢ this is the correct result 
(Lenz 1929), although without further reasoning we could have been 
confident only of an upper limit. To obtain higher order terms cluster 
corrections must be added to (2.12). This has been considered by Temkin 
(1960) with the assumed form 


f=1—Lexp[-n(r—9)] ci ae NES CE 


forr>c. The best 4 is density dependent, approaching zero (as k,) at low 
density so that (2.19) then agrees with (2.15) near the core. Only a rather 
poor upper bound is obtained for the coefficient of the c* term (see § 3.5) and 
the dependence on the number of spin states is wrong. Temkin conjectures 
that the fault lies mainly not with the assumed form (2.19) but rather with 
the Jastrow wavefunction (2.1) itself, and that generalization of the latter 
along the lines of § 2.5 is necessary. This is in contrast with the boson case, 
where Iwamoto (1958) and Aviles (1958) obtained a close upper bound on 
the first correction (occ**) to the Lenz term. In the boson case it is for this 
purpose necessary to consider not merely the next cluster term but an 
infinite subset of higher terms. 

’eturning to the general case, the following difficulty has been observed 
by Emery (1958). If the potential contains a sufficiently strong attraction 
outside the core, the solution of (2.14) does not minimize (2.13), which is, in 
fact, unbounded below. ‘To see this we suppose the modified Schrédinger 
equation 


— SV (,W)+dy=Edyh . . 2. . (2.20) 
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has a bound state }, with energy &<0, which we normalize so that 


[IvPar=1. 


Then, if fis any function satisfying the required boundary conditions, eqn. 
(2.2), (f+ By) will also satisfy these conditions where B is any constant. 
Putting this correlation function in (2.13) and using (2.20) we find for the 
energy shift per particle 


AE(f+ By) = AE(f) + (constant). B+ B2é. 


Since &< 0 it follows by letting B+ co that the energy can be made arbitrarily 
low. Ofcourse no physical significance can be attached to this result since, 
with large B, the neglected cluster corrections will be large. The result 
shows, however, that, with a potential for which (2.20) has a bound state, 
the variational method has no value. According to Emery this difficulty 
occurs for any potential for which it is possible to obtain a negative energy 
shift in (2.13), but this is not in fact true. We can rearrange (2.20) into an 
ordinary Schrédinger equation with a effective potential 


1 V2d,12 

("+ Grae) 
and the difficulty only arises when this effective potential has a bound state. 
The term depending on d, in this effective potential is fairly small, but as it is 
not of constant sign it is not possible to relate this condition exactly to the 
condition for v to have a bound state. With a realistic nuclear potential it 
is unlikely that v does have a bound state, since the existence of the deuteron 
probably depends essentially on the tensor force, which does not contribute 
in (2.12). 

However, even if the Emery difficulty does not occur, it seems likely that 
the long-range part of the correlation function, determined from (2.14) 
will be too large for the first approximation (2.12) to be adequate. This is 
certainly so if we ignore the d function since in this case « in (2.15) is the 
free space scattering length which, in the singlet state, is about — 24f. 
This may be somewhat reduced by the proper inclusion of the d functions 
but we know of no detailed investigation (see the following section, however). 


2.3. Numerical Calculations 


In this section we review the numerical calculations which have been 
made using the ideas of the previous sections of this chapter. The only 
realistic calculations which we know of using the Euler equation have been 
made by Emery (unpublished). Emery used the central part of the 
Brueckner-Gammel potential (§ 1.3) acting in even states only, and solved 
the Euler equation, i.e. (2.14) generalized to allow for a spin dependent force. 
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The results at three densities are given in table 1. The binding energy at 
equilibrium density is approximately 23 Mev per particle, corresponding to 
k,x1-6f-1. According to the work of Brueckner and Gammel (§ 3.8) the 
potential used by Emery should not give any binding energy (Brueckner 
and Gammel obtain about 20 Mev per particle from the tensor force), so 
it seems that the first approximation is very bad. In view ofthe surprisingly 
small scattering length this conclusion is unexpected, and some further 
work on this point would be useful. 


Table 1. Results obtained by Emery for the solution of the 
Kuler equation 


Binding energy per Scattering length 
particle 


20-9 Mev —1-08f 


22-8 Mev —0-85f 


22-5 Mev —0-68/ 


i 
Iwamoto and Yamado 


| : 
it—Mean Separation 


0 | 2 3 4 5 6 
CFermis) 


Showing the correlation functions in the variational calculations of Iwamoto 
and Yamado (1957) and Emery (1958). 


Other calculations have been done with restricted forms of f, in which only 
a few parameters are determined by the variation process. If the cluster 
corrections were negligible this would give an upper bound for the energy, 
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although of course it might not be a close bound. The main interest of 
such calculations is in the magnitude of the corrections. 
Iwamoto and Yamada (1957b) used the form 


Flr) =0, r<Se 
=l—exp[-—f(r—c)], r>ce. 


Their potential, fitted to low energy singlet scattering, is a hard core of 
radius c = 0-6f plus an exponential Serber attraction. The only variational 
parameter is 8, for which the optimum value is 5f-! at an equilibrium 
density with r)=1-1f, when the energy per particle is —5:2mMev. That 
the energy is too high is to be expected in a method which suppresses the 
tensor force (see § 3.8). With the quoted value of B the correlation seems of 
quite short range compared with the particle spacing. Nevertheless, the 
first cluster corrections, estimated by Iwamoto and Yamada, aresubstantial. 
The energy shift per particle in lowest order consists of 48 Mev kinetic 
energy and 77Mev potential energy. The first corrections to these are 
respectively (—12)mev and (+12)Mev, or about 20%. The numerical 
approximations used do not permit a significant estimate of the percentage 
correction to the total energy. 

In this connection the work of Dabrowski (1958 a, b, c) is of interest, 
although concerned with the particular nucleus 1°O rather than with nuclear 
matter. He used potentials, fitted to low energy singlet scattering, which 
had Yukawa-Serber attractions in addition to hard cores. For f he chose a 
gaussian form outside the core. The results are shown in table 2 for the 


‘Table 2. Showing the results of Dabrowski’s (1958 a, b, c) calculations on 
16Q, The fourth, fifth and sixth columns give the results of the 
lowest order calculation, and the last three columns give the next 
order corrections. Energies are in Mev and lengths in fermis 


A 8 ro Kinetic|Potential| Total Kinetic | Potential Total 


energy energy . 

energy | energy jenergy| | ection | correction correction 
0-2 | 0-75 | 0-86 | 984 —1183 |—199 — ae = 
OO. heei2. p15 651 —716 | —65 6 —3l —25 


two values of core radius considered. They should be compared with the 
experimental values 7)= 1-1f, H=—128mev. These should in fact be 
reasonably approximated for some intermediate core radius. The first 
cluster corrections (evaluated for c=0-6f) to the separate kinetic and 
potential energies are relatively smaller than those of Iwamoto and Yamada, 
but the percentage correction to the total energy is still large, being about 
40%. 
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In the calculations mentioned the correlation functions have the pro- 
perty f<1 everywhere. Thus any degree of positive correlation, such as 
might naturally be expected from the attractive part of the potential, is 
excluded. Emery (1958) used a form free from this restriction, namely 


f(r) =9, rs, 


f(r)= {1- Sexp[-Mr—e)]] [l+aexp(—fr)], r>c 


where «, 8 and X are variational parameters. His potential was different 
in singlet and triplet states and fitted separately to both sets of low energy 
data, despite the omission of a tensor force. He chose c=0-5f. In lowest 
order the resulting energy and density were close to the required values. 
However, because of the strong positive correlation, the first cluster 
correction was large and was estimated to contribute about 20 Mev per 
particle to the binding energy. The rapid convergence of the perturbation 
expansion with non-singular potentials (§ 1.4) suggests that the large positive 
correlation found by Emery is probably spurious (in fact with the potential © 
used it may be, as discussed in § 2.2, that the first approximation has no 
true minimum), but the calculation shows the difficulty encountered in 
obtaining a reasonable value for the energy by this method. Unfortun- 
ately there are no results available to show whether the inclusion of the 
next order term appreciably alters the position of the minimum. 

The overall situation is thus rather unsatisfactory. The work that has 
been done with assumed forms for f suggests that cluster corrections are 
essential for accurate results, and may even be so large as to vitiate the 
whole method. Certainly the labour involved increases enormously when 
it is required to calculate even the next order term in the cluster expansion. 
Finally it should be noted that for realistic results it is certainly necessary 
to generalize the trial wavefunction (2.1), at least to include tensor correla- 
tions which are required in order to obtain the contribution of the tensor 
force to the binding energy. 


2.4. The Effective Potential 


We have seen that, whilst it might be possible to treat at least the most 
important part of the correlation due to the hard core by the cluster expan- 
sion method, any positive correlation due to the attractive part of the 
potential is more troublesome. It is suggested therefore (e.g. Clark and 
Feenberg 1959) that one might use an F/ factor to remove the hard core 
(2.1) and then treat the remaining part of the potential by some other 
method, e.g. by perturbation theory. To see how this might be done we 
consider a complete set of unperturbed (fermi gas) states of the nucleus © - 
and define the set ‘l’, by j 

Y~,=f®,,. amet y', Feat pre be ensk Wee acer 
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Then, in lowest order of the cluster expansion}, we obtain 


CE eB ES =T,<¢®,| ®,) 


9 22 OA L/2M{V Ei)? + 20) 2G) IP, 
(2.22) 


where 7’, is the unperturbed kinetic energy of the state ®,. Thus, in this 
order, we can ignore the F factor provided we replace the potential v(ij) by 
an effective potential given by 


1 
vert = = (VE) + af? Sie ee et 5) A228) 


The states Y’, defined by (2.21) are in general non-orthogonal but this does 
not show up in the lowest order approximation. The effective potential 
will be non-singular and can be used in low order perturbation calculations 
(cf. § 1.4) and, in finite nuclei, in shell model calculations for example. 

The function f must be chosen to satisfy, as well as is possible, two con- 
ditions; (i) it must be such as to make the higher order cluster terms 
negligible; and (ii), it must be such that the approximations used to treat 
Vert are adequate, i.e. if perturbation theory is used then vere must be 
sufficiently ‘weak’. In practice it is far from clear, however, how one 
should apply these conditions. One possibility, which is of interest because 
it gives results which are related to those obtained by direct perturbation 
theory methods (see §3), is obtained in the following way (Clark and 
Feenberg 1959 and Moszkwoski and Scott 1960). Consider a potential v 
interacting only ins-states and separate it into two parts, a long-range and a 
short-range part, defined by 


%=0 and vg=v, r<d, Bes 
vyj=v and vg=0, r>d. \ (2.26) 
The distance d is chosen so that the solution of 
1 
— U V (d, Vf) =e > Unda f = 0, ° . . . (2.25) 


where the vs are the short-range parts of the v, (2.5), defined in an analogous 
way to vs in (2.24), areequaltolforr>d. With this/the effective potential 


is given by 


Veti =U, r>d 
; (2.26) 
=F (vpptoss, rsd. 


cy Ste ee ee a eel ee a 

+ Note that since we are not merely considering expectation values in the 
ground state we are here using the term ‘ cluster expansion ’ in a generalized 
sense. For reasonably low excited states the conditions for the first approxi- 
mation to be adequate should be approximately as before. 
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Equation (2.25) is the Euler equation for the minimization of AJA (2.12) 
when the potential is just vs. The energy shift to lowest order in vert is 
immediately seen to be just the expectation value of ) in the state ®. 
In second order there are terms in 7)? and also contributions from the part 
of vere With r<d. 

If we now make the approximation of putting /(a) =0 in (2.25), the equa- 
tion becomes the zero energy Schrédinger equation for the potential vs. 
Thus vs has been chosen so that, approximately, it produces no s wave 
scattering at zero energy. The value of d is approximately If, so the 
neglect of the higher order cluster terms is not unreasonable. 

Moszkowski and Scott (1960) use a method which is equivalent to a 
generalization of the above procedure, in which the f’s depend on the state 
of the particles concerned (see §2.5). The distance d is then a function of 
momentum defined so that vs, in the modified Schrédinger equation, gives 
no s-phase shift at any energy. Second-order perturbation theory is used ; 
the method then becomes equivalent to the Brueckner method (§3) with a 
certain approximation to the solution of the K matrix equation, and we shall 
therefore treat it further in§ 3.5. We note however that, from the viewpoint 
of this chapter, the choice of d, and hence of f, is rather arbitrary. Ideally 
one should choose d so that the sum of the neglected cluster terms and the 
corrections to the method used to treat v) is a minimum, but this is hardly 
practical. Allowing d to become a function of momentum, as above, 
certainly has the effect of weakening v; as the momentum increases (d is an 
increasing function of momentum) and hence improving the convergence 
of the perturbation expansion in 2, as pointed out by Moszkowski (1959), 
but it is not obvious that this is not at the expense of increasing the magni- 
tude of the neglected cluster terms. 


2.5. The Bethe—Goldstone Equation 


The physical ideas underlying the methods of this section and those used 
to justify the methods of §3, are not dissimilar, and there has been some 
interest in the relationship of the two methods. It has not been found 
possible to obtain the approximation of § 3 by a variational method, but we 
shall indicate here how this approximation arises from a generalization of 
(2.1) by the neglect of certain cluster terms. Some of the ideas of this 
section have been used by Leigh (1958), by Weisskopf and de Shalit (1959), 
and by Dabrowski (1958c). How far one can regard this section as a 
‘derivation’ or merely as a series of steps leading to a known result we leave 
to the readers. 

It is necessary to generalize the wavefunction (2.1) by allowing the 
internucleon correlations to depend on the states occupied by the nucleons. 
We denote the single particle (plane wave) states by ¢; (r;) and define the 
non-antisymmetrized ground state ¢ by 


Date WH) kere cB eBoy 


occupied 
states 


Theory of Nuclear Matter 235 


Then we define an antisymmetrizing operator by 


Ey gel OA eres ne, ha (2:28) 
The generalization of (2.1) is obtained by putting 
Po Aan Me Wen ee ae egal = 2929) 
with hey) wee etic Sos) & ee 2°30) 
i>j 
Si =Sii- ° ° . . ° . . . . ° (2.31) 
We thus have 3? different functions f,,. We require these to be such that 
HY = HY’ +higher order cluster terms. . . . (2.32) 
To see how this can be done we have 


H¥ =AHFS 


= {T ade [a Vif + ee(ky— ~ I) ). Why) fy | 


Vis: a 
5 OS SAEED Cy ee eee ET ORE 

aM fifa 1 | 
We are here denoting the momentum of the state 4; by k;. Each gradient 


operator acts on the argument of the f factor immediately following it. 
The j,; are then chosen to satisfy the two-body equations 


rf 1 a A 
{ MeCifi+ 5p Vf pe) - Vis ~ Qi) fy} B67) = 0 
(2.34) 
where .%,, is the operator which antisymmetrizes between particles 7 and Jj, 
and where Q;; is a projection operator which we discuss further below. 
The necessity for the f,; to depend on the states of the two particles is seen 


immediately from this equation, as pointed out by Dabrowski (1958 c). 
Equation (2.33) becomes 


Vfy. V 
HY =| AT,+4 S Ae( i) |e Al 43 S(Q'y—Velis) + we Bie sais | a. 
pea (2.25) 


If we can neglect the second bracket then we obtain, as our approximation 
to the energy 


BEA UTR SS OAT VI cart aie omen (OIRO) 
2A ij 
Returning now to the two-body eqn. (2.34) we introduce the two-body 
wavefunction 
big t) = Lig Fig (Tag) Bi (ts) $y (5) Daisy eando cri eH 


then (2.34) becomes 
AettiMby=— (3p V+ eC) ~ Qu’ Wi) du eee > (286) 
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where V acts on the relative coordinate r;; and where 


e(ij) = alk- k; [?, 


the relative kinetic energy of the Sera nucleons. We now make 
explicit use of the well-known suppression of free scattering in the interior 
of a fermi gas by the exclusion principle (Weisskopf 1950) by choosing 
Q;;' so that Q,,’ v(ij) does not scatter into occupied states. Thus we put 


(,(t) ba(9) [Qi = Lifk,? and k,? > k/? 
= 1 if either ~=i and B=j or a=) and B=1, 
=0 otherwise. 99. 4°.) « Gepeset oh ceerenmeemennrera! 


It follows that (Q,;’— 1)v only scatters into states which are already occupied 
by particles other than? andj. Hence 


(Qy’— lye} fg b= 0,0 ae eee 


so that the term containing (Q,,;’— 1) in (2.35) is of second and higher order 

in the correlation. This is also true of the remaining term in (2.35), so 

the use of (2.36) for the lowest order expression for the energy shift is 

plausible. It should be noted that these ‘cluster corrections’ are here 

wavefunctions, rather than expectation values in the unperturbed ground 

state, so the justification of their neglect is less clear here than previously. 
From (2.38) the energy shifts are given by 


Ae(ij) = (bij 09) Hig) ee ee ee (241) 

where b= %i;b:%;, <i 1s =1, re tt. meted 
Qv(ij) 

and pig = bij + au) + Ae(i) - ve (2.43) 


where Q is the projection operator onto states above the fermi surface. 
It follows from (2.43) that Ae(ij) is O(1/A) for a large system. The total 
energy shift, from (2.38), is then O(A) as required. In the infinite limit 
Ae(ij) disappears from (2.43) and the equation is then, apart from one 
aspect mentioned in the next paragraph, identical to the Bethe—Goldstone 
equation (Bethe and Goldstone 1957) which can be obtained by a partial 
summation of perturbation theory diagrams (see §3.4). The energy shifts 
Ae(ij) are (apart from the point discussed in the next paragraph) given by 
the same perturbation expansion as the diagonal matrix elements of the 
Brueckner K matrix. We shall not therefore discuss the solution of this 
equation any further here. We note however that the Q’ operator was 
chosen to some extent arbitrarily and that the choice made is a good one 
only if the solutions of (2.43) are reasonable. In fact it turns out that the 
solution to (2.43) has an unusual behaviour in the neighbourhood of the 
fermi surface (i.e. k,~k;~k,) for potentials of interest and that the neglect 
of Ae(it) in (2.43) leads to a singular expression for the total energy shift. 
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This suggests that the choice of Q,;’ may not be suitable. The problem of 
singularities is discussed further in § 4. 

There is an important difference between the theory of this section and 
that usually used in perturbation theory methods—we have here nothing 
analogous to the self-consistent energy denominators (§3.7), the e(ij) in 
(2.43) being just kinetic energies. The fact that self-consistency effects 
seem to be important in nuclei (§§3.6, 3.7) implies a serious defect in the 
method of this section and indeed of this chapter (with the exception of 
§2.4). Inorder to obtain self-consistency effects it is necessary to go beyond 
the lowest order in the cluster expansion and it seems to be necessary to 
introduce correlation functions f which are non-commuting operators (so 
that the correlation between two nucleons depends on their instantaneous 
states rather than on the states they occupy in the fermi gas). The task 
of dealing with such correlation functions appears to be very formidable. 

Finally we mention the work of Tagami (1960) who used a wavefunction 
of the form (2.30), and calculated the expectation value of H in lowest order 
of the cluster expansion (i.e. the procedure of § 2.1 generalized to allow for 
the state dependence of the f;;). He did not minimize this quantity but 
chose the f,; to have general features suggested by the Bethe—Goldstone 
equation. With a potential containing a hard core and central Serber 
attraction, chosen to fit both singlet and triplet low energy data, he obtained 
a binding energy per particle of 18-3 Mev at an equilibrium value of 7, equal 
to 0-94f. He evaluated the next order cluster corrections and found them 
to be small. With a similar potential, however, Moszkowski and Scott 
(1960), using methods described in §3.5, found very differcnt results— 
much too low energy and no equilibrium at a reasonable density. 


§ 3. BRUECKNER THEORY 
3.1. Introduction 


The qualitative ideas of Brueckner theory have been described in § 1. 
They will now be developed in a formal way by rearrangement of a complete 
perturbation expansion. The ordinary time independent perturbation 
method is not suitable here (§3.2).. The necessary time dependent formal- 
ism is outlined in§ 3.3. Then the simplest form of Brueckner theory is stated 
and corrections to it discussed. In the light of these the ‘self-consistent’ 
theory is introduced, and the work of Brueckner and Gammel described. 
Some questions, relating to ‘single particle energies’ are discussed in § 3.9 
and § 3.10, and finally some remarks made on the applicability of the re- 
arrangement method. 

Our account is based entirely on the ‘linked-diagram’ expansion. The 
simpler versions of the theory have also been obtained by rearrangement of 
ordinary perturbation theory (Rodberg 1957, Tobocman 1957) but this 
approach is dependent on the linked diagram demonstration that certain 
difficulties (§ 3.2) cancel. These same difficulties were inherent in the many 
early non-perturbative developments (Brueckner et al. 1954, Brueckner 
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1955a, Brueckner and Levinson 1955, Eden and Francis 1955, Eden 1956, 
Watson 1956, Bethe 1956, Brenig 1957, and Kummel 1957). Animportant 
current approach to many-body problems is by way of Green’s functionst 
(see for example Martin and Schwinger 1959 and Klein 1961). These 
occur incidentally in the present work (§3.10), but the infinite chain of 
coupled equations satisfied by them can be made the basis of the theory. 
The characteriste approximation then is to break off the infinite chain at 
some point, near the beginning. This is only roughly equivalent to natural 
partial summations in perturbation theory (see Martin and Kadanoff 1961). 
The method has not been applied in a practical way to the nuclear matter 
ground statet, and is not developed here. 


3.2. Perturbation Theory for Many-body Systems 

The characteristic difficulty in the quantum perturbation theory of many- 
body systems can be explained as follows. Consider a large number N of 
similar non-interacting systems. Suppose that to the hamiltonian of each 
a perturbation proportional to a parameter A is added, the systems remaining 
uncoupled. For each system the projection of the perturbed on the 
unperturbed wave function is of the form (1—cA*) for small A. For all V 
systems considered as a compound system the projection of the perturbed 
on the unperturbed state is then 

(1 —cA?)* x exp (— Ned?) x 1— Ned? + 4.N(N —1)c7A*+.... 

Even for a small perturbation the successive terms increase rapidly for 
large NV, and individual terms have no limit as Noo. At the same time 
the perturbation series for independent component systems may converge 
very well. Now in dealing with a large volume of nuclear matter we may 
(very loosely) think of it as made up of very many smaller volumes which 
interact only weakly. Then it can be expected that perturbation expan- 
sions of some quantities will contain terms which have no limit for mfinitely 
large volume. On the other hand, expansions for local quantities, such as 
energy density, should remain sensible. The problem is to extract these 
expansions from formalisms that start always with the state function of the 
whole system. 

When the Rayleigh—Schrédinger or the Brillouin—Wagner perturbation 
method is used (Brueckner 1958 makes a detailed study) terms appear in 
each order involving various powers of the total particle number. It is 
only after extensive cancellation among these that sensible results emerge. 
If only the first few terms of the expansion are required this is no great 
inconvenience, but for an analysis of the whole perturbation series some 
method of generating the sensible terms directly is desirable. Such a 
method exists already in field theory, which is concerned always with 
systems of infinitely many degrees of freedom. As applied to the present 
subject the essential ideas were set out by Salam (1953). We follow below 
the later independent approach of Goldstone (1957). 


{ The method of Coester and Kiimmel (1960) is somewhat related. 
{ See however the recent work of Puff (1961), Annals of Physics, 18, 317. 
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3.3. The Linked Diagram Series 


Let H, E and |0) be the hamiltonian, ground state energy and ground 
state of the system respectively, and let H, 7) and |0) be the corresponding 
unperturbed quantities. The perturbation we denote by V=H-—H,. 
Then the energy shift due to the perturbation is 


AE =E—T,=(0|H—H,|0)/<0|0) 


” 
Ss ae. (3.1) 
To obtain |0) we use the time dependent S matrix defined by 
S(t,t’) =exp (t Hot) exp [—iH(t—t’)]exp(—iHpt’). . . (8.2) 
Then 
Pee rt ai ih g:9) 


{0 |0) t>— 0(1—7Ze) {0 | S(0, t)| 0) 


where « is some positive number. To see this we have, in terms of the 
complete set of eigenstates |n) of H with energies H,, 


S(0,t)= > exp ((,t) |n)<n|exp (—7H,f), 


whence the right-hand side of (3.3) is 


In) exp [i(Hn—Tp)t] (n]0) 
ibe Bh Iisa Ee ees Fal id gy A ewe Bie | 
o-sato SO n)Pexp Ea Tol S77 


In the limit the terms in (3.4) with the lowest value of H, (i.e. H,=H)=£) 

dominate, so, provided A is such that the unperturbed ground state is non- 

degenerate (a ‘ closed-shell’) and provided (0|0)40, the result follows. 
From (3.2) the S matrix satisfies 


£ S(t,t')=-i VO S60) Pee R Ne ee (3:5) 
where 
Vit)=exp (Hot) V.exp(—18,f)--. ... . . . (3.6) 


Operators transformed as in (3.6) are said to be in the interaction representa- 
tion. With the initial condition S(t’t’) = 1, the solution of (3.5) by iteration 
gives 


t t ty 
S(t,t')=1+(-4) { a Va) + (—i)? } a i “dl, V(t) Va). 
= LP Paty, TW) Vila) Vl) OD 
n : t 


In (3.7) all time integrations are over the same interval, the symbol 7 
denotes that the V(t) must be rearranged in order from left to right of 
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decreasing t and the 1/n! factor compensates for the fact that the same 
contribution is counted 7! times in the integral. 
Then, from (3.1), (3.3) and (3.7) 


oly fos af dt dt, T(V(0) V(t). “10? 
AK= Lt 
t>— 0(1—ie) wpy ee tT (Vine eee Vt,)) 10) 
: (3.8) 


We now use the notation of second quantization. The unperturbed 
hamiltonian is written 


H,=- far dt(e) Vaud), <> 5-4 eee 


57 
and the perturbation 


Ve 1 [ae dr’ st(r) bt(r’) o(r—r’) d(r’) b(n). rr oak} 


The latter form is for a local interaction, but the generalization is trivial. 
The field operators are defined by 


=F u(r) A, 
k 


(3.11) 
= > 4.* (ra 
k 
The u, are a complete set of one particle wave functions, with kinetic 
energy eigenvalues k?/2M, with periodic boundary conditions in a cube 
of volume Q. For spinless particles they are just plane waves 


uw, = 2 SFaxp Gk Flow cen ae see Seem ene 


with each component of k an integral multiple of 27Q-'3. For particles 
with spin and isospin we consider the label k to specify a spin state as well as 
a momentum, and the uw, to contain also spin wave functions. The states 
u, and operators y then have spinor indices which we normally suppress. 
In (3.10) for example v becomes a matrix in spin space and there are sup- 
pressed summations over the spinor indices of the field operators. The 
creation and destruction operators, a,* and a,, satisfy the anticommutation 
relations 

[ays dye] = [4 7, eT], = 9, \ 


Srmaeed bs) 
[a Ue] = Sic } ( ) 


A destruction operator acting on the vacuum gives zero. The many-body 
unperturbed ground state is 


[0)=Ile+| vacuum ),"\— Joa aces oe ae) 


where as before the label m indicates a state below the fermi surface. 
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Finally we introduce the interaction representation field operators 
(r,t) =exp (tH ot) f(r) exp (—iH yt) = ¥ exp (—ith?/2M) Uj Us 
k 


We*(rr, t) =exp (iHyt) W(r) exp (—iHyf) = > exp (ith?/2M) u,*a,+. (3.15) 


The final expressions here follow from (3.11) and (3.13) with (3.9) in the 
form 
Hy=)> (k?/2M) a,+ a,. 
k 


The interaction V(t) can now be written 
Vij= 1 far dr’ ht(r, t)bt(r’,t’) o(r—r’) b(r, t) (r,t). (3.16) 


With (3.16) substituted in (3.9), evaluation proceeds by means of Wick’s 
theorem (see, for example, Hamilton 1959). To state the result it is 
convenient to regard the time arguments of the ¢ in (3.16) as infinitesimally 
earlier than those of the 4+. Then the symbol 7’ in (3.8) can be regarded as 
specifying also the order in which field operators are to be taken within a 
single V(t). It is also convenient to extend the definition of 7’ to imply a 
change of sign when an odd permutation of fermion operators is required in 
passing from the written to the temporal order. Then the expectation 
value of a 7'-product of field operators is a sum of terms corresponding to all 
possible ways of associating in pairs the % operators with the + operators. 
If there are n operators of ack kind there are (2n)! ways of pairing: The 
contribution of each such pairing is the ee 

CO[L(ble, #) P(e’, #'))] 0) C0 | LHe", tyre", ")) 0)... (3.17) 
of the expectation values of individual ie multiplied by a factor (—1) 
if the written order of operators in (3.17) differs by an odd permutation from 
the original written order. The various contributions to the numerator 
and demoninator in (3.8) for a given value of n are conveniently represented 
by diagrams constructed as follows. We draw n vertical dotted lines 
labelled 1 to n to represent the n operators V(t,)...V(t,), and, for the 
numerator, a further dotted line labelled 0 to represent V(0). We refer 
to these lines as ‘interactions’. Then we join the ends of these interactions 
by directed solid lines such that one line enters and one line leaves each end 
of every direction. Each of these lines represents the pairing of two field 
operators. There are (2n)! different ways of doing this for the denominator 
and (2n+2)! ways for the numerator. The contribution of any diagram 
obtained in this way is 


0 
=| dt,.. ty {Ody dy DY LO, Vi cbge cs XasY pxta) tes) 
’ t . 
for the denominator and 
0 
=| Aly Aly | Ag Ay Dey. AY nL Oo VosiX teen) (3.18) 


P.M.S. i" 
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for the numerator, where x, y are associated with the upper and lower ends 
of the interaction respectively. The integrand J is the product of the 
following factors: 


(1) A factor 
(20)? v (Ky, —Yn) e e e © e ° ° (3.19) 


for each interaction n¥ 0, and a factor (2)! v(x, — yg) for the 0 interaction. 


(2) A factor 
Q(x, 2") = O[T(W(x),w*(x"))[0). . . . (8.20) 


for a directed line from an upper end w’(= x’,t’) to an upper end 2, and 
similar factors 


Gy, y'), —@Q x,y’), —Pyjz') . >. 2 (3:21) 


for lines from lower to lower ends, etc. Actually the minus signs here always 
cancel for the complete diagram. 


(3) The directed lines form a number of closed loops, for each of which 
there is a factor (—1). 


The lowest order numerator diagrams are 


0; of! ewan eahics (eae 


In the next order there are 4! = 24 diagrams of which a few are 


0} I ret sexs arate) 


Cy ee am 
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Wee 
SOS 


ISA she ty tn (08) 
| l cep a Mewes 27) 


| Meet iss 2 (3,28) 


The diagrams (3.24) and (3.25) illustrate that the general numerator 
diagram consists of a part connected to the interaction line 0, and one or 
more separate parts. Now the sum of numerator diagrams containing 
such disconnected parts cancels exactly against the denominator in (3.8). 
For consider the general numerator diagram with interactions 0, 1,.... 
Denote by « the part containing the interaction 0 and n(«) others; let I, be 
the integrand which would be appropriate in the absence of other parts. 
Similarly let I, be the integrand for the remaining part B, considered 
separately as if it were a denominator diagram of order n(8). We need 
consider explicitly only those diagrams in which « involves the particular 
interactions 1, 2, 3,... n(x), multiplying by a factor (N !/n(«) ! (8) !) to allow 
for the other possibilities. Then the general numerator contribution has 
the form 


N!} a) 
na) nip)! Jk GE [axe dy, dx, dy, dt,...1, (Xo, Yo U1 Yr b13 +++) 
rae ae 
Summing over all possible diagrams («, 8) the total numerator is 


(—1) ) ron) 
ae Tyrant | Modte-- 1, | dx, dy. dy. or 


={ >! — [exo Lh] SAT [do dad. . ,\. 


The second factor is just the expansion of the denominator, so we have the 
stated result. The quotient (3.8) is then the sum of numerator diagrams 
all of which are linked to the 0 interaction. Symbolically 


AE=<0|[> = i ; % Ab eed TV (OV Ves V by) One mat a28) 


where L means that only linked diagrams are included. 
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The propagator G® is explicitly 
G9(x,0') = (0 | T(x), P*(a’)) 9) 
= {0K -v’) > -at-t) > bn u(x’) exp[ —i(t—t')k?/217] 


Ikl >ky Ikl <kf 
= a {at > -at-t) > jexp [ik(x —x’) —i(t—t’)k?/2 M1] 
Q Iki> kg Ik! <kp 
(3.30) 
or, in the limit Q> ~, 
(an )= (2n)-2{ -1)| dk — 6(t' -o| ak 
k>kg k<ky 

x exp [tk(x—x’)—a(f—U’)k2/2M)]. . >. . . (3.31) 


With respect to spin indices G® is the unit matrix. 

Provided that G® falls off sufficiently rapidly for large differences in its 
arguments, it is clear that for a linked diagram the integrations over 
X,Y, ¢, etc. converge to a definite value for Qc (we assume here that v 
has a finite range). The integration over x, and y, then gives a single factor 
Q, so that the energy shift per unit volume approaches a constant. Thus 
every term in the linked diagram expansion for AZ behaves reasonably. 
A diagram with n disconnected parts would give an integral proportional to 
Q”. It is the explicit cancellation of such terms at the outset that con- 
stitutes the solution of the problem stated in § 3.2. 

From a given linked diagram distinct but equal contributions can be 
obtained (a) by permuting the labels 1... of interactions other than 0, 
and (b) by interchanging the roles of opposite ends (wy) of one or more of 
these interactions. Only one member of such a class is ever considered 
explicitly ; the remainder are allowed for simply by omitting (a) the factor 
n! in (3.18) and (0) the factor (4) for each interaction other than 0. 

Using (3.30) the space integrations in (3.18) can be replaced by sums over 
all values of labels k, one attached to each directed line. An interaction 
such as 


Ki Fay T 
| 
| 
| 


k, —_—<—_1____<____k . 
is then represented in the summand by a factor 


(k,'k,’ |v [ky ky) = | Uk dy tye,* Uy. *0(K—Y) Uy Uy, (3.32) 


This must contain a momentum conservation factor 


Stk,/-+ke’), (ky+k,)* 
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Because of this many diagrams that can be drawn do not in fact contribute. 
For example the diagram (3.26) 


oe 
— 


kZ 


2. 


could contribute only when k, and k,’ correspond to equal and opposite 
momentum. On the other hand, it follows from (3.30) that it could con- 
tribute only with k, and k,’ on opposite sides of the fermi surface. There is 
therefore no contribution. In general we do not draw or mention such 
diagrams further. 

When (3.30) is used the time integrations also can be performed. A 
diagram with (n+1) interactions gives n! distinct contributions, one for 
each possible order of the times ¢,...¢,. It is convenient to represent these 
parts separately by time-ordered diagrams, as distinct from the original 
Feynman diagrams. In the new diagrams the order of the interactions 
becomes significant; those with earlier time argument standing further 
to the right. The associated integrals run over all positive values of the 
intervals between successive time arguments. For example the Feynman 
diagram 


breaks up into the four time-ordered diagrams 


—— —— 
— ee 
and two more which vanish by momentum conservation. It is clear from 
(3.30) that a line directed from right to left must be associated with a state I 
above the fermi surface and a line from left to right with a state m below. 
Let > (/?/2M/) be the sum of kinetic energies of lines crossing from right to 
left an interval between successive interactions and let > (m?/2M/) be the 
corresponding sum for lines crossing the other way. Then on such 
interval contributes to the time integration a factor 
. a 

TS OMe ea e883 

[__ dep SPeM-Se P= ssa O88 
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This time integrated form is usefully compared with the result of ordinary 
time independent perturbation theory. The Brillouin—Wigner series 
for the energy shift is 


Q Q Q 
ee dil ee ee era ENG, | ee a St) 
A= O1V+V a7 oH, A ee I?) 


(3.34) 


where, as before, 7’, is the unperturbed ground state energy. The pro- 
jection operator Q excludes the unperturbed ground state. When expanded 
in terms of unperturbed intermediate states the series (3.34) can be 
represented by diagrams like the time-ordered diagrams introduced above ; 
the forward and backward lines represent the particles and holes in the 
virtual excited states. To obtain from this the linked diagram expansion it 
is necessary to make the following changes (which compensate one another 
exactly) : 


(1) The quantity AZ is omitted in energy denominators. 


(2) Terms represented by diagrams with parts not connected to the 
extreme left interaction are omitted. 


(3) The matrix elements of V are computed ignoring the presence of any 
particles or holes other than those actually scattered, created, or 
destroyed, by the particular interaction involved. 


As a result of (3) there are non-vanishing contributions from intermediate 
states in which several particles appear to occupy the same state above the 
fermi surface, or several holes the same state below}. Thus the lines in the 
diagrams for the linked cluster expansion strictly do not denote occupied 
states or holes, but only operator pairings in the application of Wick’s 
theorem. 

Implying both (2) and (3) in the subscript L we write symbolically 


Bah ant fae 
= (O1V + eso 
{0 car joy} 


The. existence of such a linked diagram expansion was conjectured by 
Brueckner (1955 a,b). It was established by Goldstone (1957) essentially 
asabove. The introduction of time is not of course essential in a stationary 
state problem. It is avoided in the much more complicated work of 
Hugenholtz (1957a). Several other papers also claim to establish linked 
diagram expansions in a time independent way (e.g. Riesenfeld and Watson 
1956, Coester 1958), 


| This point is much emphasized by Katz (1960 a, b). 
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Of course the method of this section reproduces the low order perturba- 
tion formulae used in§1. In first and second order the only contributions 
are 


(3.35) 


| : Oe Pee EG BS 
eke? 


These lead at once to (1.19) and (1.33), the second member of each pair 
being the exchange contribution. 


3.4. The Low Density Limit 

As has been observed in § 1 individual terms in perturbation theory are 
not sensible with hard core potentials, but reasonable results can perhaps be 
obtained by summing an infinite subset. The subset proposed by 
Brueckner can be given a formal basis (following Hugenholtz 1957b) in 
the low density limit. Now nuclear matter is not a low density system, the 
overall range of the forces being comparable with the interparticle spacing. 
However, in so far as the inclusion of higher order terms is dictated mainly — 
by the very strong short-range part of the potential, in the selection of 
such terms low density considerations are relevant. 

According to Hugenholtz the dominant time-ordered diagrams at low 
density are those with the least number of backward, or hole, lines. He 
argues that each such line is associated with an integration over momenta 
below the fermi surface and so introduces at low density a factor k,’. 
Actually momentum conservation often cuts down the number of such 
integrations which are independent, and clearly only independent inte- 
grations below k, should be counted. The argument would then be 
sufficient if the remaining integrations were over a finite domain and if the 
integrand were bounded independently of k,. The first condition may be 
effectively secured by the fall-off at high momenta of the interaction matrix 
elements—always supposing for formal arguments that the hard core is 
very large but not actually singular. However, the second condition does 
not hold; the integrand is unbounded because energy denominates can be 
arbitrarily small. Moreover, in at least one case the formal sum of the 
series has a low density behaviour different from that attributed by this 
argument to the individual terms (§ 3.5). Nevertheless, for want of better 
arguments, we proceed on the assumption that a first approximation need 
include only those terms with the minimum of independent hole lines. 
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Every diagram has at least two independent hole lines. Those with only 
two are of the following types, with any number of interaction lines: 


Plater Cattapy fees see et EY 
SATEEN ESS 
Meise 


Note again that diagrams differing from these only by the turning over of a 
number of interaction lines, such as 


or 


are automatically allowed for by the modified rules'of calculation on page 244. 
With one interaction line only we obtain on the two diagrams the first 
order energy shift 


4+ > (mm’ atest Ar i ceathor ae ate 
It is readily confirmed that all higher order ‘ladder’ diagrams (3.36) and 


(3.37) are included by replacing the matrix elements of v by those of 
K(m?+m'?/2M) where K is determined by 


v [IV | K(&) 


-} (Pa) PM . 
=0+ 2 (py M te 
or, in more compact operator notation, 
K(é)=v+v—2 _ K(6). 
(6)=v+e o—w,*() chatty aye re ian: Veraeetsy 


Here #, is the two-body kinetic energy operator, and Q is a projection 
operator onto states above the fermi surface. Solution of this equation by 


iteration : 
é)=v40-,2 PEN ow 
AC Tet PORTS peor nN Aer: 


Ur noes Mor GAO) 
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provides explicitly the contributions of individual higher order diagrams. 
Equation (3.39), together with 


AEZ=34 > (mm’'|[K (“Sr ) [mm’—m'm), . . (3.41) 
ant 2M 

defines the simplest form of Brueckner theory. It is obviously correct to 

second order in perturbation theory and contains all those terms which 

appear to be important for a strong core of sufficiently short range. In the 

following section we discuss methods for solving (3.40) and applications 

thereof to nuclear matter. 


3.5. The K-matrix Equation 
The expression (3.41) for the energy shift involves directly only the 
diagonal and exchange matrix elements: 


m? +m’ m2+m’2 


«{mm’ |i “S) |mm’) and (mm’'|K S) |m’m). 


These are ‘on the energy shell’, in that the argument & of K(&) equals the 
total kinetic energy of both initial and final states. We consider first an 
approximation in which such matrix elements are expressed directly in 
terms of scattering phase shifts without reference to the potential. 

We have to solve 


K(€)=v0+L(&) K(é) Sica mr Sect Bx4e) 
where 
axa 
=) so Ee Bae (eae 
W6)= > 3-11 RM oe) 
For comparison, the free space reaction matrix K° is given by 
KG H=V vb Gy ayy te es © (3,44) 
where 
| kk’) kk’ | 3 45 
(O)=P % geek DM ee) 


Here the intermediate states are not restricted to being below the fermi 
surface. For the interesting values of & the denominator can then vanish, 
and the principal part of the singular integral is taken, as indicated by the 
Pin (3.45). Now K can be expressed in terms of K°, for a formal solution of 
(3.42) is 

Kerk eke i-th) Kop eae pa (8.46) 


If the free space scattering K° is weak the second term can be ignored. For 
spinless particles then 


(mm | K((m?-+m2)/2M1) | mm’ )=Q> MED (21-41) tan 8, (k) 
(3.47) 
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where 5,(k) is the scattering phase shift for the /th partial wave and k is the 
relative momentum 

k=}|m—m’|. 
Similarly 


(mm | K((m?+m’?)/2M)| m'm) = Q-4 it > (21+1) (—1) tan §,(k). 
1 
(3.48) 


For nucleons it is convenient to work, not with the independent spin 
states of two particles, but with the singlet-singlet, singlet—triplet, triplet— 
singlet, and triplet—triplet combinations. In the summation (3.41) these 
take the appropriate statistical weights 1, 3, 3, and 9 respectively. For 
each combination, assuming central forces only, formulae (3.47) and (3.48) 
hold with appropriate phase shifts, except that (3.48) changes sign when 
the overall spin state is odd under exchange. The energy shift is then 


2a 2 
eS ool es 2 13 31(]; 
AE o>. 1z,2, 1+1)(3 tan 6%(k) +3 tan 634(k)) 


: > (214+ 1) (tan 6,(k) + 9 tan spr) RN Sag hls. 
Zodd 
With 

4Q 
A =4 = 
> sai|am 


in the limit Q-> oo, and the integration over (m+ mm’) carried out, the 
energy shift per particle is 


AE 4k2 (rk dk Bo DAuN 
A aM}, i (=) } 


x1 > (21+ 1)(3tan $43 tand%)+ > (21+ 1)(tan $+ 9 tan 65%) 
leven l odd 
(3.50) 


For nuclear matter at normal density the phase shifts needed are those 
up to a laboratory energy of about 160Mev, the most heavily weighted 
being around 50 Mev. 

Formula (3.50) and a generalized version permitting tensor forces were 
used in the early work of Brueckner et al. (1954) and Brueckner (1954). The 
phase shifts were computed from a potential determined in part by meson 
theory and with a phenomenological hard core. Reasonable values were 
obtained for the energy and the mean spacing when the tensor force was 
omitted, but its inclusion resulted in both quantities becoming much too 
small, 

In such calculations there occurs the difficulty that, because of the bound 
deuteron state, the triplet s phase shift approaches 7 rather than 0 at zero 
energy, and passes through (7/2) in the range of integration of (3.50). The 
integrand is then singular, and Brueckner et al. took the principal value of 
the integral. However, in such circumstances the formula is just not 
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applicable. The derivation implied weak scattering and so small phase 
shifts. Indeed we are entitled to make the substitution 


UNE > Oat ie rata Torr ecole weno 35D. 1) 


which at least removes the singularity if the formula is used for large phase 
shifts. Now it happens that (in the absence of bound states) such a formula 
then gives the exact answer for AZ in the trivial ‘many-body problem with 
one-body forces’ (Bell 1959, Katz 1960a, b). Moreover, the second virial 
coefficient for a non-degenerate quantum gas with two-body forces is also 
linear in the phase shifts (Beth and Uhlenbeck 1936), with again extra 
bound state terms. So attempts have been made to provide in the present 
case also some real reason for preferring 6 to tané (Fukuda and Newton 
1956, de Witt 1956). That these are futile is indicated for example by the 
case, now to be considered, of the pure hard core interaction. There either 
6 or tané gives correctly the term of first order in the core radius, but 
neither predicts the presence of the second-order term. 

For the pure hard core the s-wave phases are simply kc, where c is the 
core radius. So from (3.50) to first order in ¢ 


(=) -(AF): Bo tnnnecben ta 5) 


a formula due to Lenz (1929). Consider now further terms in a systematic 
expansion in powers of (k,c). The expansion of tan6é for s-waves gives a 
contribution of order (4,c), as do higher partial waves. Diagrams other 
than the ladders (3.36) and (3.37) are of third and higher order in perturba- 
tion theory, and so also of at least this order in K°, or c, in a rearranged 
series}. The only contribution of order c? is then from the second term in 
(3.46) which is adequately approximated by 


RCO Ee LOI ty Se ete 23 (3°83) 


To compute this, off-diagonal matrix elements of K° are required. It is 
readily shown (by methods such as those used later in this section) that to 
lowest order in c 


; 1 [4ac 
{m,m, | K(&) | my! m,’) = Q-1 5 (m, + mg), (my' + m,’) (Fr) j 
(3.54) 


Using this in (3.53) and (3.41) and performing intermediate state integra- 
tions (Levinger et al. 1960) yields the correction to (3.52) 


ke \a0 6 

soe 11=—2In2) (k,c)?: rom SiO) 

ot ( ~) 3 aaa ) (ky 3 
; 

+ On dimensional grounds they then contribute as (ksc)®, which in terms of ky 
is of order ky? relative to the leading term. On. the other hand, such diagrams 
have one more hole line at least, and on this ground might have been thought 
to be of relative order k/?. This is the point mentioned at the beginning of § 3.4. 
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This result was obtained in essentially this way by Martin and de Dominicis 
(1957) and Galitskii (1958), and in different ways by Huang and Yang 
(1957) and Yang and Lee (1957). Martin and de Dominicis gave also an 
approximate value for the next term. Including this, we have 


AE _ kf 66.63 (hc) + 0°32 (kc)? + 0-78 (kc)? + ...} (3.56) 
tary 2 ve ; oT oi ; 
For nuclear matter, with k,=1-5f-t and c=0-4f7, 
AB/A =40{0-35+0-1040:14+ ... } Mev. 
Clearly the core radius is not small enough for convincing convergence. 
The terms given add to about 24 Mev, about the same as the mean unper- 
turbed kinetic energy. 

If the core radius c were sufficiently small and the remainder of the 
potential sufficiently weak, the latter could be allowed for simply by adding 
its expectation value in the unperturbed state to (3.56). This approxi- 
mation was used for nuclear matter by Gomes et al. (1958). Their potential, 
fitted to low energy singlet scattering, was a hard core plus square well 
attraction. The whole potential, including the core, had a Serber exchange 
character. The coefficient of the third order term in (3.56) with no p-wave 
contributions and no contributions from diagrams other than the basic 
ladder diagrams (3.36) and (3.37), was found to be 0-26 rather than 0-78. 
The equilibrium spacing came out very close to the empirical value and the 
binding energy per particle was about 7Mev. ‘This result was thought 
reasonable in view of the omission of the tensor force. 

The validity of this approximation was examined by Levinger ef al. 
(1960a, b). Let vg denote the hard core and 2 the remainder of the 
potential, so that v=vs+v (this corresponds to taking d=c in (2.24)). 
If Kg is the solution of (3.42) with vs alone, K can be developed in powers of 
K,andv,. As far as second-order terms: 


K=K,+y4+Lm+KgLu4+uLbKs. . . . «+ (3.57) 


The first two terms are those retained by Gomes et al. (1958). To estimate 
the corrections Ay, is replaced by the first approximation (3.54). The 
third term is just the second-order energy due to a regular potential and 
might be expected to be small from the work of §1.4. However, the pre- 
sence of the hard core requires v; to be stronger than would otherwise be 
necessary to fit low energy scattering. Asa result the second order term in 
v1 is quite large, contributing about (—20) Mev per particle to the energy. 
The remaining terms in (3.57) express the fact that the hard core reduces 
the wave function in its vicinity and so partially suppresses the effect of x). 
Levinger et al. find from these terms an energy of about (+20) Mev per 
particle. There is considerable numerical cancellation among the correc- 
tion terms, but they are individually so large that the convergence of the 
method is very doubtful. 

Levinger ef al. used for 7, conventional forms in which the strongest part is 
also the innermost, and so most affected by the proximity of the core. 
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With the square wells of Gomes et al. the correction should be smaller: 
The numerical examination of this particular point by Moszkowski and 
Scott (1960) is vitiated by their omission to vary the effective mass self 
consistently with density. 

The method of Moszkowski and Scott (1960) exploits the substantial 
cancellation between the effects of the hard core and the innermost part 
of the attractive potential. By allowing vs to include in addition to the 
hard core some part of the attraction, out say to a definite radius (see 2.24)) 
it is possible to weaken both K,° and v and so reduce the second-order 
terms. One choice of d is such that vs gives no free space scattering at zero 
energy. Not much appears to be gained by letting d vary with momen- 
tum and angular momentum. Moszkowski and Scott estimate the second- 
order corrections} for a Serber potential, with 0-4f hard core and an 
exponential attractive well, fitted to low energy singlet data. The last 
two terms in (3.57), the cross terms, are the largest ; they contribute about 
2Mev per particle at normal density. The other terms are less by an 
order of magnitude. Thus the method is much more acceptable than that 
of Gomes et al. With the potential mentioned no equilibrium is found near 
normal density. However, by omitting the substantial d-wave and small 
g-wave contributions (i.e. by retaining interaction in s-states only) reason- 
able equilibrium density and energy are found{. It is argued that this 
actually simulates the effect of the tensor force. But clearly it is desirable 
to extend the method to handle the tensor force properly, and also the large 
self-consistent field effect found by Brueckner and Gammel (§3.7). This 
has now been done (Moszkowski and Scott, private communication) §. 
Essentially, the results of Brueckner and Gammel are reproduced with 
much less labour. Details are not yet available. 

Finally we outline the direct numerical solution of the K matrix equation 
by partial wave analysis. The spin states are dealt with in the singlet— 
singlet, singlet—triplet, etc. combinations; a separate equation holds for 
each, involving the relevant part of the potential. Because of conservation 
of momentum, the total momentum p of the two particles is a fixed para- 
meter. We write 


ip+k’, Ip—k’ | K(6)|4p+k, p—ky= Ck’ |K(E,p)[k) (3.58) 


where <k’ | K(é, p)| kk) = <k’[o|k) 
elk” f(k", p)(k’|K(8, p)| k) 
Re ae ee Grd Se 3009) 
k pit tae ol Pl US as oat Mal —k’|2 


In these estimates some allowance is made for self-consistent modification 
(see § 3.6 and § 3.7) of the energy denominators, but this should not alter the 
order of magnitude. 

+ In these calculations the third order self-consistent field effect of § 3.6 is 
included, and plays an important role in locating the equilibrium. 
§ See Bull. Amer. phys. Soc., 6, 149 (1961). 
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and f expresses the exclusion principle effect 
f(k’,p)=1 if [4ptk’|)k, Rperemer vr cad 
=0 otherwise. 
Clearly 


(k’ | K(8, p)|k) = [arexp ik’ ror), (7)... (3.61) 
where #,(r), dependent of course also on & and p, is determined by 


(ir) =exp (ik .r) + | dr’ G(r, r’)u(r’)b,(r’) . . . (8.62) 
with 


2 2)-1 
G(r, isa p) exp [¢k. (r— ‘i — ra = ar . (3.63) 


Following Brueckner and Gammel (1958) we work directly with the integral 
eqn. (3.62). Actually the first general investigation of the K-matrix, made 
by Bethe and Goldstone (1957) used an integro-differential equation 
obtained by operating on (3.62) with 


Lic 2 
(a ae t 6 — i) . 

This gives the inhomogeneous Bethe—Goldstone equation 

(e- pay o) (r)= (s- arte (ik.r) + [ar g(r, r’) v(7’) ob, (r’) 

4M M or?]'* 4M M ; 
(3.64) 

with g(r, r')= > (k, p) exp (tk. (Fe) Or ee eee (3.65) 

k 


The equation becomes homogeneous on the energy shell, & = (p?/4.M + k?/1) 
Equation (3.62) differs from the free space scattering equation only in 
that G contains the exclusion factor f. This has two important effects. 
The first is that the energy denominator is never zero (for values of & 
needed in (3.41)), and the scattered wave therefore decreases more rapidly 
thanr-!. This is important when ¢, is considered as a Jastrow correlation 
factor (§2.1); it is of no special importance here, as is clear from (3.61). 
The second effect is that for p40 the Green function @ is anisotropic. Ina 
partial wave analysis the various waves are then coupled even for a central 
potential. ‘The general coupled wave analysis is set out by Werner (1959). 
Brueckner and Gammel remove the difficulty by approximating f by its 

average f over directions of p. The quantity f is unity if 
| cos @ |< (kp)—1 (4p? +k?-k?) . 2. 2 2 6 © ©(3.66) 


where @ is the angle between k and p, and zero otherwise. Averaging 
over angles we obtain 
f(p,k)=0, tp +k<k?e 
=1, |[}p—k|>k;, 


=(kp) (4p? +k?—k7?) otherwise. ©. . (3.67) 
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When % is then resolved into partial waves (characterized by quantum 
numbers L, S, J, J,, 73) the only remaining coupling is that due to the 
tensor force between: triplet states with L=J+1. The radical functions 
for these states are determined by pairs of integral equations, those for 
other states satisfying separate single equations. The details of this, 
given by Brueckner and Gammel, raise no question of principle. The 
only other point we consider is the special treatment demanded by the 
hard core. 

To illustrate this we consider for simplicity only the singlet s-wave. 
The s-wave part of exp (ik.r) is 


afe du exp (tkrp) = —— a (meat ements ee. a0) 
From (3.62) the appropriate integral equation is 
+ | “Grjr'yo(r')dlr)dr’ se. (8.69) 
with : 
Gnrj=arr | feveau ao 
4M M 
=, (alee Mee ae ae emia. bar tort) 
Ser tiatar 


Now with r<c, the potential is infinite. From (3.69) ¢ must then be zero, 
and the otherwise indeterminate combination vé must be chosen to bring 
this about. To simplify the problem note now that the inside of the hard 
core has no physical significance, and can be removed if we wish without 
changing physical results. An extreme case of such a ‘hollow’ hard core 
is a 6-function of infinite strength atr=c. 'The contribution of this to v¢d is 
of the form Aé(r—c), where A is a constant to be determined (it represents 
the product of the infinite strength of the potential and the zero wave 
function). Then (3.69) becomes 


$(r) = ¢9(r) +AG(r, ¢) + | Td Car woNbr)  < (3.71) 
which together with 
$(c) =0=¢9(c) +AG(6, ¢) + | dr’ Glo, yur \h ©), (3.72) 


gives 


+ | * ar’ {(r.7)- Oj r')d(r’). . « (3-78) 


The hard core is thus eliminated at the expense of complicating the Green 
function and the inhomogeneous term. This equation and corresponding 
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ones for other waves are suitable for computation. The A-matrix contri- 
butions are finally obtained as integrals of the type 


[arguinrrer gy 
=go(c) + | * dr bor) (7) $(0) 


ae Ry ye se G(c,r)\ 
ane | au {u(r bl) Te 3) o(r) d(r). (3.74) 
Now this approach gives an accurate solution for what may be called 
a ‘thin hollow hard core’. It is also an approximate solution for a ‘solid’ 
hard core, but there corrections are needed because (dv) does not vanish for 
r<c. There are indications (Bethe and Goldstone 1957, Brueckner and 
Gammel 1958, Appendix B) that such corrections are negligible in nuclear 
matter. Any appreciable difference in results, for differing assumptions 
about the inside of the core, would indicate inadequacy of the approximate 
theory being discussed. 

The method just described was applied by Brueckner and Gammel 
(1958) using the potential described in §1.3. The Green functions were 
constructed and the integral equations solved numerically, with rather 
coarse lattices. Actually the work was performed merely as the first step 
in a much more elaborate programme (§3.8). The main point of interest 
is that the binding energy per particle from this first approximation 
exceeds by about 20 Mev (at normal density) the final result. So with this 
potential, the most realistic used in such calculations, the simple version 
of the theory so far developed is rather poor. 


3.6. Third-order Corrections 


We consider now corrections to the Brueckner ladder approximation. 
All first and second-order perturbation diagrams are already included, 
but most of the third-order diagrams are omitted. With singular potentials 
the perturbation theory terms are not immediately useful, and we go over 
to a modified expansion in K rather than v. The low order terms are most 
quickly generated simply by solving (3.42) for v in powers of K, 

v=K(€)—K(€)L(€) K(é)+... Ret nahgict ac'eg eee 
and substituting in the perturbation series. The first corrections to the 
theory of the last section are then of third order in K, and obtained merely 
by using matrix elements of K rather than v in the third-order perturbation 
diagrams. Note that any reasonable value of & is permissible in (3.75) ; 
the é& dependence of K is clearly a second-order effect and is ignored in the 
crude estimates to be considered, The use of K rather than v is roughly 
equivalent to considering not individual third-order perturbation diagrams 
but whole families generated from them by replacing interaction lines by 
ladders. 'The cavalier treatment of the & dependence of K prevents such a 
correspondence being precise. 


Theory of Nuclear Matter 257 


A particular class of third-order diagrams is obtained by the insertion, 
into one or more lines of second-order diagrams, of self-energy parts. 


Examples are 
° ! 
(in 2 een (oas) 
Cn 


ane 


This class we set aside for the present. The remaining third-order diagrams 


are 
—— ‘particle—particle’ scattering (3.77) 
ae > ‘particle-hole’ scattering Gee S18) 


ee, “hole—hole’ scattering ¢ Wei seat 37D) 


where, to avoid for the moment distinction between their ends, we have 
used interaction lines of zero length. The first of these is actually included 
already in the ladder sequence. ‘The others are new. 

The ‘ particle-hole scattering’ is also known as the ‘ three- -body cluster’ 
term; this is because it has three hole lines and may be thought of as a 
collision of three particles from the fermi sea. It was estimated by 
Brueckner *(1955b), with extensive revision by Bethe (1956). The 
component diagram 


I m I m f 
Pa een Maa pat Siete a Jegs@a, (5.50) 
ear 
gives a contribution 
_ <m' m|K|V1)¢m'l | K | m’ ae m). (3.81) 
(m2? +m’? — 2 —1'2)(m? +m"? — —1"*)/4M? 
U 


P.M.S. 
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The factor (4) normally associated with the first interaction is omitted ; 
this is compensated by a factor 2 arising from the alternative time ordering 
of the other interactions. Distinct diagrams are obtained by letting a 


line to the top of an interaction go instead to the bottom and vice versa, 
for example 


Such variations can be allowed for by the substitution in (3.81) 
K— K(1—P'P’P’) 


where the P are exchange operators in space, spin, and isobaric spin. We 
must not consider also interchanges of lines going away from an interaction ; 
this would involve counting separately diagrams differing only by inter- 
changing the roles of opposite ends of interactions. As regards the second 
and third interactions this is already allowed for in the rules of calculation. 
As regards the first interaction it is equivalent (in the Feynman diagram) 
to interchanging the remaining interactions, also allowed for already. 

Brueckner and Bethe approximate the matrix elements of K by those 
of a conventional Serber potential fitted to low energy singlet scattering. 
The approximation cannot be good; it does not for example reproduce the 
repulsive effect of the hard core which dominates for large momentum 
transfer. However, it should give an order of magnitude. With w just a 
function of internucleon distance then 


KS ye (GER UE Ss sist ges tae ea aes 
K(1—P*P°P")=4w(1+P)(1—P°P”). . . . . « (8.83) 


The spin summation implied in (3.81) involves the trace over the spin 
states of three particles : 
Trace 
(1 —PyoP2)(1 — Pog’Po3")(1 —Pgy°Psy7)=12. . . (3.84) 


This is calculated by using the expressions 


P,.°=3(1+6,.0,), etc., 


4] 


and reducing the expression to one linear in the spin operators of each 
particle, whose traces are zero. Brueckner (1955 a, b) quotes a spin factor 
24. Bethe (1956) obtains the factor 12, but thought that he had omitted 
some terms (see his footnote); however, they are just exchange terms 
already allowed for. With the notation 


l=m-—q, =m‘+q, I"=m’"+q 
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and converting sums to integrals, the total contribution to the energy shift 
per particle is 
9M? 
Tir ls 
4 OD + (14 + m= m])}{w(g) + v(m" — m’ |)}{e(q) + (|g +m" — m))} 
(P?+q.m—q.m’)(7?+q.m—q.m’) 


(AE/A)pp = jaa dmdm'dm’ 


(3.85) 
where w(q) is the Fourier transform 


w(q) = fer exp (—7q.r)w(r). 


Bethe finds}, with an approximate evaluation of the integral, 
(AH/A)pp + — 0-3 Mev 
at normal density. Brueckner (1958) has a similar result. 

With the particular approximation used to evaluate the integral, the 
correction is, at low density, proportional to k;, and so more important 
relative to the leading term in the energy shift (ock,3). Bethe suggests 
that this indicates the breakdown of the Brueckner method and the 
formation of real clusters. However, the approximations used in the 
integration were suitable only for high density (Brueckner 1955b). The 
result, not very reliable at normal density, therefore permits no conclusion 
about low density. According to Hugenholtz (§3.4) the cluster diagram 
would be relatively less important at very low density because it has an 
extra hole line. But this would not mean the improved convergence of the 
method as a whole, for clustering and indeed condensation must occur with 
realistic forces. These phenomena escape the essentially perturbative 
approach. 

Bethe (1956) has made a second estimate, to some extent complementary 
with the first. In this the hard core alone is considered. In the short 
range limit (from (3.54)) 

w(q)=4re/M. ae eee eee TSO) 
In this limit the assumption of Serber exchange character is unimportant 
(P’= 1) and we can substitute (3.86) in (3.85). Simplifying the integral 
by ignoring the m? in the denominators and integrating over all ¢ greater 
than k, we get the rough estimate 

(AE/A)pp © 0-021 (kc)? (k,?/2I1) 

~ + 0-2 Mev 
with c=0-4f, ky=1-5ft. 

Finally we quote a value for the cluster correction at normal density, 
due to Brueckner and Gammel (1947) : 

(AH/A)pp > +0-3 Mev. 
The context suggests that this has been obtained from a K matrix resulting 


+ Actually he quotes a somewhat smaller value due to the use of an effective 


mass (§ 3.7). 
U2 
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from numerical solution of the appropriate equation with a realistic 


nuclear potential. 
The hole-hole scattering diagram (3.79) has direct and exchange com- 


ponents 


Together these give an energy shift 
(AB)nn=4> (my mpl K[hils) tC my m,’)¢m,'m,'|K|m,m,— m,m, ) 
(my'? + mg’? — 1,2 — 19?) (1my” + m4” — 1,? — 1,")/4.M? 

(3.87) 
This could be treated much as (3.78). However, the only estimate in the 
literature proceeds in another way. Consider the expression 


m2 Mo" wor me oie m2 Mo 
(mimil (5 2 ae )lm, m,’)<m,’m, KS ag Siz) mim. mam,), 


(m,? +m" — my? —m,")/2M 
(3.88) 
where the singularity is excluded in the summation. Although apparently 
of second order in K, this is actually of third order, for, if the & dependence 
of K is ignored the summation vanishes identically because of the anti- 
symmetry between the m’s and m’’s. To second order one readily sees 
that 
Br Sica eee Q 
K(8)-K(6)=K (zo - ee Z)K . . B89 
where #, is the two-particle kinetic energy operator. Using this (3.88) 
is seen to be equivalent in third order to (3.87). Moszkowski and Sessler 
(1960) evaluate (3.88) for the A matrix resulting from the work of 
Moszkowski and Scott. They find an energy shift per particle 
(AH/A hh ~ 0-5 Mev. 
It is of some interest that the hole—hole scattering diagram can be included 
along with the particle—particle scattering in a sequence of diagrams known 
as ‘generalized ladders’. These are of the form 
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where the interactions (other than that at zero time) may occur in any 

relative order). We discuss later (§ 4.3) the summation of all such diagrams 

due to Prange (1960). According to Iwamoto (1959, 1960) an approximate 

Jeena is effected by replacing the original K matrix by K” determined 
rom 


K"(@)=v+v so Kae Ee aera (S00) 


Here P projects below the fermi surface. Equivalently 


—P 


K"(8)=K(6)+K(6)z—~ K"(6) 
oer 0 
=K(6)—K(8) gp K(6)+ KO) gy KE) gy K(é) 
0 
BP e 


The term apparently of second order here is just the third-order term of the 
preceding paragraph. Indeed Moszkowski and Sessler started from 
Iwamato’s equation. Fortunately they did not proceed to the explicitly 
third-order term, which is fictitious. 

It remains to consider the insertion of self-energy parts into second-order 
diagrams, as in (3.76). The insertion may be into any of the four original 
lines. Constructing the corresponding energy shift in perturbation theory, 
but with matrix elements of K rather than v, yields 


(AB) se, 
; E\hle) (ile | K| mm, 
cee ee A mae? (VC) + V(a)— Voom) — Vemma) 
(3.91) 
where 
Vik)= 3 (km|K[km—mk) . 2... (3.92) 


m 


The quantity V is a generalization of the effective one-body potential 
of Hartree-Fock theory, which is obtained with v rather than K in (3.92). 
To the order of accuracy considered the term (3.91) may be incorporated in 
the original formula (3.41) by the substitution 


m?/2M > m?/2M + V(m), (3.93) 
2/2M + I?/2M + V(l), hy hp ene aa 


in (3.39) and (3.41). This means in effect that a particle propagates always 
in a mean field, due to other particles, which modifies its energy. This is 
what we have meant by occasional references to ‘self-consistent field 


+ Like Moszkowski and Sessler we assume here that singular integrals are to 
be defined with principal values. Iwamoto’s intention is not clear and his 
derivation suggests that the point is one of some delicacy. 
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effects’. If the matrix elements of K are replaced by those of a regular 
central potential the energy shift (3.91) is quite small. According to 
Martin and de Dominicis (1957) it is no more important than other third- 
order effects. However, with more realistic potentials the result is rather 
larger. Moszkowski and Scott (1960), with a central Serber potential with 
hard core, find a shift of about 5 Mev per particle at normal density. Witha 
still more realistic potential (that of Brueckner and Gammel 1958) their 
estimate greatly increases (private communication). It is then in sub- 
stantial agreement with that obtained by Brueckner and Gammel, as 
quoted at the end of §3.5. The increase is attributed mainly to the 
inclusion of the tensor force. 

The general situation then is this. In so far as corrections have been 
estimated, only the self-consistent field effect appears to be substantial. 
This certainly cannot be omitted and may be too large to be adequately 
treated as a small correction. The latter view is taken by Brueckner and 
Gammel and motivates their approach in which self-consistent field effects 
are built in at the start. To this quite elaborate formalism we now turn. 


3.7. The Self-consistent Field 


The systematic introduction of the self-consistent field depends, in the 
language of field theory, on the elimination of self-energy parts. We work 
here entirely with Feynman (as distinct from time-ordered) diagrams. A 
self-energy part in one such as 


SEINE): neath 


which becomes disconnected from the rest of the diagram when two solid 
lines are broken. It can be regarded as the result of inserting additional 
pieces in a single line. The ‘skeleton’ of any diagram is obtained by 
removing every portion that can be regarded as such a self-energy insertion, 
with the proviso that the special interaction at t= 0 must not be so removed. 
The family of diagrams with a given skeleton is obtained by making all 
possible insertions independently into all solid lines of the skeleton By 
associating with s a line origins anti , 
es g i ; such a line not the original quantity G(x, x’), but a new 

quantity G(x, x’), obtained by summing over all possible insertions, the 
contribution of the > fami ‘diagrams is giv indi 
8 aeiecre : e ae fi amily of diagrams is given at once. To indicate 
7 rather than G" is implied, we draw skeleton diagrams with heavy 
lines thus 


_ 
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Self-energy parts are called ‘ proper’ if they cannot be divided into two or 
more parts connected only by asingle line. The diagram (3.94) is improper, 
being a succession of three proper parts 


The contribution to G(x, x’) from proper insertions only we write in the 
form 


[eter der are, 0" (x" | iM |x") G@%x",a').  . . (3.95) 


Here the factors G® arise from the terminal lines and (—iM) is a sum over 
all possible proper insertions. It is usually convenient to regard the 
variables x, x’, etc., as generalized matrix indices and to suppress them as in 
_ matrix algebra. Then (3.95) is written simply 


G° (—1M) G®. neta dene pt ass yee 90) 
With this notation the complete propagator, the result of making any 
number of proper insertions, is clearly 
G=G—G91MG+GiMGiMG—.... 
This is the iterative solution of either of the integral equations 
oie. \ eee a Nie (9.07) 


It is useful also to eliminate ‘ladder’ parts, i.e. those of the type 


The sum of such ladders with any number of interaction lines we denote 
by a single heavy interaction line, 


| fT 
i] 
1 


ee 


To obtain the appropriate factor note that the basic diagram 


xl! 
| 
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makes a contribution which can be written 


[ery de” ay” G(x, oa) Gy, y")<x" y" | —w Lary" > (2™, x’) Gy”, y’) 


(3.98) 
where 
(ay |v |a’y!) = v(x — y) B(x — Yo) (Xo! — Yo’) 8(9 — Xo) 5(K — x’) B(y—y’) 
(3.99) 


where x, is the time component of x, etc. The next contribution 


: xl yt ! 
| ' 
I ! 
PS SE ee 
4 y!! yl! y! 


can be obtained by making in (3.98) the replacement 


«xy | —iv |a’y’) + [at dy" dx” dy” 
x {ay | —tv|x"y") G(x", x")Gy",y”) Kx" y” | —w|a'y’) 
or in an obvious matrix notation 
—iw—>(—w)GG (i). 


In this notation the sum of arbitrarily long ladders is obtained with the 
replacement 


—iv—> —iK’' 
K'=v+0( —1G@)v + o( —iGG) v( —14G@) v 
ato rays 


The matrix K’ is thus the solution of either of the integral equations 


K' =v—-wG@@é K’, 
K' =0-iK’ GQ. } Sates! 
The quantity K’ is related to the Brueckner reaction matrix, but the relation 
is complicated by the fact that the latter was defined for ladders in time- 
ordered rather than Feynman diagrams. To spell out the connection note 
first that in all orders (ay | K’ |x’y’) contains the factor 8(ay — yp) 8(ao’ — Yo’) 3 
so there is really only one initial and one final time argument. Define 
retarded and advanced parts of the propagator G by 


(3.101) 


where 
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Define a new quantity K by 
K=v0-wG+G+ K 
=a OA GON EN vaals Ci umanl sot. (aeL02) 
This will have a momentum representation of the form 


(ay | K |2x'y’) =8(a9— Yo) 5(aq’ — yp’) 
> <kike | K(%o, %9’)| ky’ ky’) exp (—ik, . x —iky.y +ik,/ x! + tk,’ .y’). 
Finally introduce the transformed operator K(@) by _ 


K(é)= | dtexp(—iét) K(0,t). . . . . (3.103) 


Then it is readily seen that in the approximation G@=G@°, K(&) is identical 
with the operator previously so labelled (3.39). Note from (3.100), (3.102) 
that 
K'=K—-ikG-G-K’ 
age hae Oh wyatt, cays (3: 104) 


To first order in K, K'=K. This approximation is much used in the 
following. 
In analogy with (3.35) the first diagrams of the new theory are 


and the corresponding energy shift is 


AE=4 | dx dy dx’ dy’ (xy | K" |x'y! (Glx', «) Gy’, y) —Ga',y) Gly’, x)} 
(3.105). 


with x)=Y)=0 and integration over only negative values of x’, yy’. The 
corrections to this will not include the large self-consistent field terms, 
already incorporated. The other third-order corrections will have some- 
what modified values, but should remain small. 

The eqn. (3.105) does not become useful until some practical approxi- 
mations for K’ and Gare defined. Here we derive the equations used in the 
work of Brueckner and Gammel (1958). This derivation is much indebted 
to the work of Klein and Prange (1958) and Thouless (1959 a, b). 

The principal approximation is to calculate M only to first order in K. 
Then G is computed from M by (3.97), and K from G by (3.102). It is also 
expedient to make the approximation that G, like G°, is translationally 
invariant in time—i.e. a function only of the difference to the two time 
arguments. This property does not hold accurately because the time 
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origin occurs explicitly in the formalism, as the upper limit of all time 
integrations. This assumed time-translation symmetry of G is communi- 
cated automatically to K and M. ; 

In this scheme JM is obtained from the lowest order diagrams 


with the approximation K’=K. ‘The result is 


ce | fe’ y= fay dy’ Gey y ){Geu LK [x'y' d= Cou LK Lye" Yh 
(3.106) 


The momentum representation of this and similar quantities (diagonal 
by space-translation symmetry) is of the form 


{x | M| x’) => exp [ik. (x—x’)] M;, (xo, %9’). 
k 
Like K, M is a retarded operator, i.e. zero for <2). Because of the time 
translation symmetry (in the present approximation scheme) 
M ,(t,t’)=M,, (0, t’ —t), 


so M,, is entirely determined by the transform 


Mys)=[ dt M,,(0,t)exp(—iét). . . . (3.107) 


—c 


In momentum space (3.97) becomes 
GAé, t) = Gt, t’) ~i far dt” Gt, t") M,(t",t") G,(t",). (3.108) 


For momenta below the fermi surface, and for t’ = 0, it can be confirmed by 
substitution that the solution is 


Gin(ts 0)= — 4(—t) exp (—?e,,f) lat de fe =e (3.109) 
where 
m2 
Cn = aM + is i a - . . : ; ; (3.110) 


For ¢’ < 0 we use the translated expression 
G,,(t,) = — Ot’ —t) exp [—ie,,(t-t’)]. . . . (3.1094) 


Now this is not a solution of (3.108). It solves an equation (Klein and 
Prange 1958, 'Thouless 1959 a, b) in which the restriction t” <0 is replaced 
by t’<?t’. This mutilation can be defended on two grounds: (a) It would 
in effect occur automatically if K(t,t’) and so M(t,t') contained a factor 
d(t—t’), or equivalently if K(&) were independent of &. The error in 
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pretending that this is so is of second order in K. (b) Because one inter- 
action occurs always at zero time, the time argument ¢=0 is of particular 
importance. For this reason the error in (3.105) due to the use of (3.110) is 
actually of fifth order in K. 

For momenta above the fermi surface G;° is, like M » aretarded quantity. 
As a result the time integrations in (3.108) do not extend up to the origin, 
and G has time translation symmetry without further approximation. 
It is also a retarded operator, and so fully determined by the transform 


aye) | een CaCnOwe, ele. SO) 
From (3.108) 
G(E) = G P(E) +1G(6) M(é)}4 
={@—PP2M—-Miey-. ©. | (8.112) 


The result (3.109) has a simple interpretation. For propagation below 
the fermi surface, in this approximation, the effect of the medium is to add 
to the kinetic energy a momentum dependent potential with matrix 
elements J,,(e,,). The dependence of 1, on & in (3.112) prevents the 
inversion of the transform from giving a simple exponential for Gj. As 
usual this & dependence is of second order in K, so that G{(t, t’) is approxi- 
mately exponential. 

Next we calculate K from (3.102). Taking the transform we obtain 

TACN =O WING VI G)s a Gee hol ae (S119) 
where 


0 
(WF L(E)] 1’) =| dtoxp (—16t)G,(0,t)G/(0,#), (3.114) 
and L has no matrix elements involving momenta below the fermi surface. 
This integral can of course be expressed as a convolution of the transforms 


(3.112) of the G’s. | However, if we approximate G, by an exponential the 
integral is trivial. Suppose that e, can be chosen so that 


A,(t) = G,(0, t) — exp (ve) a ere tre (OsL ko) 


is small. Ignoring terms quadratic in A, 


Ci’ | L(é) |W) = [ dt exp (—i6&t) {exp [i(e, + et] + A,(t) exp (et) 
ee + A,(t) exp (tef)}. 
The terms linear in A vanish if 
0={€—e,—P/2M — ME —%)}*-{@-a—ey" 
+{6 —e,—12/2M —M(e —e)} *-{ 6 — ey — ej}. 


Allowing e, to depend on & and 1’ we have the solution 


e(l’, €)=P/2M+M(E—e(,€)), . - - + (3-116) 
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and likewise with J, l’ interchanged. That e,should depend on the additional 
arguments indicates simply that the best exponential fit depends on where 
the propagator is to be used. The final equation for K is 


oe ae [Wy <u’ x 
K(é)=0+0 3 gr Hayne Rae en bs 


Using (3.109) to rewrite (3.106) as, 


M,(€) => (mk|K(E+e,)[mk—km), . . (3.118) 


the circle of eqns (3.117), (3.116), (3.118) and (3.110) is complete. These 
equations are equivalent to those in Appendix A of Brueckner and Gammel 
(1958). They appear to have arrived at them (actually at the less compact 
version in their main text) in a largely intuitive way. 

A simpler set of equations is obtained by using once more the slow 
variation of M(@) with (&). Equation (3.116), differs from (3.110) only 
in the argument of M. Apart from an error of second order we can take for 
all k 

€, =k?/2M + M,, (e;,) =k?/2M + V,.. ee rer tek 


Then e, depends on J only, not onl’ and &. Only the argument & =e, is 
now needed in M,,(&). and the matrix elements of K required in (3.118) are 
on the energy shell. For k<k, the K matrix equation now involves a 
vanishing energy denominator, and the principal part of the singular integral 
is specified. The quantity V,, is a generalization of the Hartree-Fock 
potential. We call this approximation that of the ‘effective field’®. It 
differs from the previous approximation by a term in the energy which is 
fourth order. 
A further simplification once popular is to approximate V;, by a quadratic 
form 
V,=a+dk?. Cee dd Se aera 
Then 
€, =a+bk? + k?/2M 


Sb FASE on eal se eee eT 


where M* is the ‘effective mass’ (Weisskopf 1957). In the K matrix 
equation only differences of e,,’s occur, and the only effect of the modified 
propagation is to substitute J/* for the original mass. For nuclear matter 
this effective mass approximation is poor (Brueckner 1955 a, Martin and 
de Dominicis 1957 a). Typically V;, turns out to have the form shown in 
fig. 2. 

Important transitions go from below k; to well above, where V,, is 
negligible. ‘The effective mass approximation, chosen for example to 
give the slope correctly at k;, greatly overestimates the increase due to V of 
the corresponding energy denominators. 
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Finally, from (3.105) and (3.110), the energy shift is 
AEB=%3 ¥ (mm'|K(e,,+¢,,)[ mm’—m’m) (3.122) 


m 


The last equation is as in the Hartree-Fock theory; the total energy 
shift is one half the sum of ‘level shifts’ of individual occupied states. 


3.8. The Brueckner—Gammel Calculation 


Attempts to apply the self-consistent theory have been made mainly by 
Brueckner and his collaborators. Here we pass over early attempts using 
simplified forces and the phase shift approximation to the K matrix 
(Brueckner and Wada 1956, and references therein). We come at once to 
the work of Brueckner and Gammel (1958), the most ambitious attempt on 
the nuclear matter problem yet published}. This uses a quite realistic 
potential in conjunction with a simplied version of the theory contained in 
equations (3.117), (3.116), (3.118), (3.110) and (3.124). To recapitulate, 
the energy shift is given by 


Nis (Gna aa ye et. (8.125) 

where ‘ 
Cm = m?/2M + > {mm'| K(e,,+e,)| mm’—m’m), . (3.126) 
Fees AEDES) (3.127) 


rr 6 —e(l', &) —e(I, &) 
e(l’, &)=2/2M + ¥ Um] K(6 +n —e (I, &))[Im—ml). (3.128) 


To simplify these equations Brueckner and Gammel replaced{e,,, — e,(l, &)} 
in (3.128) by some ‘suitable’ average value (—A). The quantity e,(/’, &) 
defined by (3.128) then becomes independent of J’. We use the new 
notation 

el’, &) > &(&, 1) 
so that 
é(€,1)=2?/2M +> <Im|[K(@—A)[Im—ml). . . . (3.129) 


The matrix K(&) is then required for & equal to all values of (€,, + €,,— A), 
i.e. from (2e,—A) to (2e,—A) where é) and e, are the values of e,, with 
m=0andm=k, respectively. Variation of A from zero to (¢;—¢€))—about 
70 Mev—increases the equilibrium spacing by about 4% and the energy 
per particle by about 1mev. Brueckner and Gammel believe A= (e,—e) is 


+ However, calculations have been announced by Martin and Puff (1960) 
and by Moszkowski and Scott (private communication) which use the same 
potential and get similar results, apparently with less labour. 
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about right and quote final results with this value. Lacking their insight, 
we cannot see why twice this value (for example) should not be equally 
reasonable, and we would attribute a corresponding uncertainty to the 
result. It should be said that this approximation, i.e. replacing (3.128) by 
(3.129), does not seem so essential now as when the work was begun. At 
that time they had only an elaborate series of equations involving an infinite 
hierarchy of different K matrices. It was later observed (Brueckner and 
Gammel 1958, Appendix A) that the whole hierarchy could be expressed in 
terms of the one matrix K(&) involving the extra parameter &. 

The numerical solution of the equations is by the obvious iterative pro- 
cess. At first the energies are given the unperturbed values 


Cm=m/2M, &é,l)=?/2M. 


The AK matrix equation is then that already considered in §3.5, and is 
solved numerically as described at the end of that section. New approxi- 
mations to e,, and é(&, 1) are then obtained from (3.126) and (3.129), and a 
first approximation (with no self-consistency effects) to the energy shift is 
obtained from (3.125). The process is restarted with the new values of e,, 
and é(é, 1). On this second time round the modified energy denominators 
would be a new cause, in addition to the exclusion principle, for coupling 
between partial waves. Again this is avoided by spherical averaging. 
For intermediate states of momentum ($p + k) the propagator is 


f(p, KS — 28, /(4p?+ +k. p)]—28, VV (dp?-+k2—k.p)]}2. (3.180) 
For small momenta 
;(¢\P 
é @ Rak ) 


is linear in (}p + k)? and the k. p terms in (3.130) cancel. This is again so at 
high momenta, where é(é, k) is essentially just the kinetic energy k?/2M. 
To make some allowance for k . p in the intermediate region it is treated as a 
perturbation. ‘To second order in k.p we then need the average over 
(k.p/kp) of 
J(p, k) J(p, kK)(k.p)?e le, V (hk? +49") ] 
S-2WE, JCP) (F-2IEV ETP” 

where e(6, 1) is the second derivative of @(&, 1) with respect to 1. With "3 
from (3.67), and 

J (p+ k) (ks p)*=3(ko)tpA ey 30S a (oa) 


and reabsorbing the correction into the denominator, we have the averaged 
Green function 
G(r, r’) 
I(k, p) exp (ik. (r—r’) 
m= 


a a ei 
i o—a 6. [(u+in'+ + inf) | -6| 6 J (e+ gp Dh, f 
4 V3 if ; 4p V3 ‘pf 
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This is the form (resolved into partial waves) used by Brueckner and 
Gammel, except that instead of f(k,p) in the denominator they have f?; 
Brueckner (1958) hasf#?. When the calculation was repeated at Cambridge 
(Clement et al. 1959) it was found that the distinction between f? and [32 was 
negligible. Presumably then f gives nearly the same results, and indeed 
perhaps this term can be omitted. 

About five iterations were sufficient tu obtain constant results. Very 

coarse meshes were used (for example e,, was considered at only five values 
of m between 0 and k,). This could be reduced to about forty-five minutes 
by using still cruder meshes and certain integration formulae. The major 
item in this numerical work was the construction of the partial wave 
Green functions. 

Their results are quoted in table 3 for two potentials B and C. Both fit 
low energy data, have hard cores, spin dependent Yukawa central forces, 
and tensor forces; C has also a spin-orbit coupling, which improves the fit 
to high energy scattering data. In the preferred case, C, the binding energy 
and equilibrium spacing are closely equal to the desired values. In view of 
the many approximations in the theory these results might be thought 
adequate. Brueckner and Gammel observe that a crude allowance for 
coulomb repulsion in real heavy nuclei improves the agreement. As 
observed in §1.1 an equally crude allowance for surface tension gives a 
larger effect in the opposite direction. 


Table 3. Properties of nuclear matter as calculated by Brueckner and 
Gammel 


Potential B | Potential C 


Binding energy (Mev) 18-5 15-2 
Equilibrium spacing 7 ( f') 0-95 1-02 
Compressibility (Mev) 167 172 

M* at fermi surface 0-72M 0-73.M 
Symmetry energy (Mev) 32 22 


Theoretical values are given by Brueckner and Gammel also for com- 
pressibility, symmetry energy, and ‘effective mass at the fermi surface ’. 
The latter, defined by 


(M)-1 = 9 (ilar?) er [aw —2 Seo ga tcko) 


has no direct physical significance (§3.9). There is no clear empirical 
result with which to compare the predicted compressibilities ($1.7). The 
empirical symmetry energy is intermediate between the two quoted values. 
The difference of these is attributed mainly to the inclusion in C of odd 
state central interaction. The calculation of symmetry energy allows 
only for the fact that the summation in (3.122) is below different fermi 
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surfaces for neutrons and protons. It ignores the dependence of the K 
matrix itself on the neutron—proton ratio. A similar approximation in the 
discussion of separation energy (§ 3.9) gives a quite wrong result. 


3.9. Separation and Rearrangement Energies 


In the Hartree-Fock approximation there is a celebrated theorem which 
states that the fictitious single particle energy of the topmost occupied 
level has a real physical significance in a large system. It is then the least 
energy required to remove a particle from the system—the ‘separation 
energy’. The question arises as to whether the quantity e, of Brueckner 
theory has a similar significance. The answer is no, as we will see. We 
come later (§ 3.10) to the generalization of the above theorem by Hugenholtz 
and Van Hove, but this does not refer to e;. The difference between e, and 
the separation energy, say e,’, been called ‘rearrangement energy ’. 

In developing the perturbation series for the energy of the A particle 
system we assumed closed shells, i.e. all states occupied up to the fermi 
surface. Suppose now that one of the topmost states is unoccupied, and 
apply the same theory. The result is tentatively identified as the ground 
state of the perturbed (A—1) particle system. The unperturbed state 
is now degenerate with many others in general (with different directions 
for the hole momentum or different spins), but because of conservation laws 
this is no complicationt. The series for the energy of the (A —1) system 
differs from that of the A system merely in that | (particle) summations 
extend over an extra state and m (hole) summations over one less. The 
difference, the separation energy, can be represented by breaking one or 
more lines of the original diagram to indicate that they refer to the special 
state and must not be summed over. In a large system it is in general 
necessary to break only one line of a diagram at a time, for each such 


breakage removes a factor Q associated with a summation. Thus the 
original diagram 


gives rise to the separation energy diagrams 


RCE esse 
ec a es 


} We have of course momentum and isobaric spin conservation; parity con- 
servation is sufficient to show that states of different ordinary spin do not mix. 


ee GHEes 
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However if two lines refer necessarily to the same state they must be broken 
simultaneously. Thus 


(3.135) es (3.136) 


We will later (§3.10) see a construction of the separation energy that does 
not involve these disconnected diagrams, but which does not give an 
explicit series. 

The series just described can be constructed essentially by differentiating 
the original energy series with respect to k,. The only difference is that 
spins and momentum directions are then automatically averaged over, 
but this should not matter in alargesystem. Thusif H(k,) is the sum of the 
original energy series 


,__. OB(k,) Ok, 


Cuorgecone on 
_ (1) 2m) 1 fake) 
i Cees me pe  @8T) 


We could in the same way have defined a quantity e,” as the energy gain on 
adding a particle to the lowest unfilled shell. By similar reasoning we 
would have arrived again at the value (3.137), so that in a large system 
e,=e,". Note that this result has been established only for perturbation 
expansions and may fail when these are inapplicable as in the theory of 
superconductivity (§ 4.4). 

When the system is at equilibrium density the energy per particle is 
stationary as a function of ky, i.e. 


O (Bxe al a Be) <0 
Fel Gaon Wee Ae 


Then (3.137) is equivalent to 
e¢ =H/A, ° ° ° ° . . e (3.138) 


P.M.S. 
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i.e. the separation energy is the same as the energy per particle. In the 
work of Brueckner and Gammel, at equilibrium density, e,= — 28 Mev and 
E/A = —15evso the rearrangement energy is 12 Mev. 

In the Brueckner approximation 


2 
B=> > +4 ¥ (mm’ | K(e,, + é,)| mm’ —m’m ) 


so, from (3.137) 


in 2 
es = ge +> (km |K(e,+e,,)| Kym — mk; ) 
“s ak, @ 


+4 ps <mm' | 0A Fea TS aetenee Re 


=epth S (mm IS ae, Ent mm’—m'm). — (3.139) 


The last term in (3.139), which is the rearrangement energy, would be zero 
in the Hartree-Fock approximation where K reduces tov. More generally 
if we write the K -matrix equation in the form 


K=v+v(Q/e) K 


we see that an explicit dependence on k, arises both from the exclusion 
principle operator Q and from the energies e,, and e, in the denominator e. 
Retaining each of these effects to lowest order in K gives, for the change in K 
due to a change in k,, 


580 1 
k= K 2K + KOS (2) K. te ea siiae) 


Using this in (3.139) we obtain second and third-order contributions 
analogous to diagrams (3.134) and (3.136), but with self-consistent 
propagators (§ 3.7) so that (3.135) does not arise. These correction terms 
were estimated by Brueckner ef al. (1960) using the K matrix of Brueckner 
and Gammel (1958). They found 


/ 
Cy —Ce= 10 Mev 


in reasonable agreement with H/A—e, as we expect. It should be noted 
that any discrepancy here has nothing to do with the basic approximation 
of Brueckner theory, within which (3.138) and (8.139) are identities at 
equilibrium ; it reflects only further approximations made in evaluating the 
two formulae (numerical, spherical averaging, the expansion (3.140), etc.). 
Various other estimates have been made of the terms arising from (3.140) 
(Thouless 1958, Squires 1958, Mittelstaedt 1958), but these have used 
crude approximations, to the K matrix for example, and are not quantita- 
tively meaningful. 

For the separation energy we have seen that the term of first order in K 
is subject to substantial correction. This raises the question as to whether 
the first-order treatment of the operator M in Brueckner theory is adequate. 
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This allows for such insertions in particle and hole lines (of time-ordered 


perturbation diagrams) as 


| | | | I 
! } | | | | 


| | | t ! I < 
O teens 
| ] Tae ' et 
—>—_1_>— 
ee bee 
It does not allow for 


ee eet tee etc., (3.141) 


Cc, 


nor 


\ ! etc. (3.142) 


These last diagrams are clearly similar to the rearrangement energy 
diagrams discussed above. But they refer now to the propagation of 
virtual particles and holes and not to the energies of real particles and holes. 
The ‘energies’ of the virtual states are not well defined, for the propagators 
are not simple exponential functions. We saw in §§3.7 and 3.8 that only 
after a series of approximations could such energies be assigned. By 
further approximations of a similar type new virtual ‘energies’ can be 
defined which allow roughly for the additional insertions (Brueckner and 
Goldman 1960). Brueckner et al. (1960) find that the original energies are 
displaced by only a few Mev (5 at the fermi surface, rather than the 12 Mev 
for the removal of a real particle). This is substantially smaller than the 
first-order self-energy effect (e,—k,?/2M x —70Mev), which when intro- 
duced into K-matrix energy denominators increased the real mean energy 
per particle in nuclear matter by 20mMev. Brueckner e¢ al. find that the 
additional ‘rearrangement’ effect is only about 1-5 Mev. 


3.10. The Optical Model Potential 


Here again we consider the addition (or removal) of a particle to (or from) 
the closed shell system. Now we remove the restriction that the particle 
is added just above the fermi surface, or removed from just below. The 


x2 
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general case presents the difficulty that the unperturbed state is degenerate, 
or nearly so, with many states in which one or more other particles are 
excited at the expense of the energy of the extra particle. The precise 
pattern of degeneracy shifts rapidly with choice of the additional plane wave 
state and with the size Q of the box. The energy shift due to the interaction 
fluctuates correspondingly and there is no definite limit for large Q for a 
state of prescribed momentum. It is then useful to study, instead of 
individual levels, the time dependent propagation of the extra particle or 
hole. 

Denoting by |0) the ground state of the system without the extra particle 
then 

yi(r)]0) 

is a (properly antisymmetrized) state in which an extra particle has just been 
created at the point r’. The subsequent evolution of the system is given by 


exp [—1H(t—’)] #*(r')|0). 


The probability amplitude for finding the extra particle at the point r at 
time t, with the remainder of the system in its ground state, is 


Gtir,t; r’, )=O(t—-')O|p(r,t)prt(r’, ’)[0). . . (3.143) 
Here @ is the usual step function, and 


w(r,¢=exp (tHt)b(r) exp (—iH?), 
wt (r, t) =exp (tHt)s*(r) exp (—iHt), \ 


are Heisenberg field operators. For convenience we are here taking the 
energy of the state [0)to be zero. In the same way the propagation of a 
hole is described by 


G~(r,t; r’,t)= —O(t’ —t)(O|pt(r’,t’) P(r, 60). . (3.145) 


The conventional ‘one-body propagator’ is the combination 
G(r,t; r’,t’)=Gt+G- 


= O(T(p(r,t)pr(r’,t’)]0). 2 2. (3.146) 


Let |N) be the complete set of energy eigenstates for the system with 
one extra particle, and | M) the corresponding set for one hole. Then 


G*(r,t; 0,0)=6(t) 2 (0| ¥(r)| N)<N| ¥+(0)|0) exp (—i2y 0), 
G>(F5 £5 0,0)= —O(—t) © CO] *(0)]M)<M f(r) 0) exp (Ey! | 


(3.144) 


(3.147) 
where Hy, and H,, are the energies of |N) and |M) respectively. These 
have the spatial fourier transforms 
Gy*(0) = 00) ¥ Ja, | N)<N |a,+ 0) exp (—iLt), 
N 


G,-(t) = —(-t) > <0 | a,* |M)<M |a, [0) exp (iByat). (3.148) 
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Introducing a small damping factor «> 0 the quantity 
G,*(t) exp (— ef) +G,-(t)exp(et) . . . (3.149) 
does not differ appreciably from G,(t) for t<e—. It has the Fourier 
transform in time 
€) oe (0 [ee INN Ie? |» » <0 bee eal OD 


(3.150) 
When ¢ 0, for a finite system G,(&) has a series of poles, the location of 
which give energies of the (A +1) particle system with momentum +k. 
If ¢ is kept large compared with the level spacing the discrete structure is 
averagedaway. This happens ineffect when in calculations the limit Q—> 00 
(sums + integrals) is taken before e+0. A similar averaging process occurs 
in the theory of the optical potential in finite nuclei (see, for example, 
Brownf 1959) to eliminate rapid fluctuations associated with individual 
compound nucleus states. The resultant propagator still contains statis- 
tical information about levels of the (4+1) nuclei. Specifically, in the 
integral representation 


pete) B,(é’) 
= ee ee Fe, ae hell 
eke) ifae ever ER SR ATE Boy 
the positive functions A and B are weighted level densities 
€)=>5(E—EHy) Ola, |N)CN[a,+|0), } 
N 


(3.152) 
= 5 5(6 — Hyg) 0 [m+ [M)<M | 4, [0), 


where the ‘6’-functions are supposed of width sufficient to cover many 
levels. 

The propagator G can be represented by a linked diagram expansion of 
the type already used. We again need the relation (3.3) between real and 
model states 

_  S(0, — )|0) 
Eis Ana (0, — c)10) 
and the similar formula 
oj= —K01So.0) 
{0 | S(co, 0)|0) 


The Heisenberg operators are given in terms of interaction representation 
operators by 


(r,t) = S(t, 0) d(r, t) S (é, 9), 
woe t)= S(t, 0) bt (r, t) S (, 0). } Sees etal o3) 
Thus ae o 
trea! <0 | S( S( oo, t)b(r, Poke 
Giirjésr,t)= ane aaa 
(3.154) 


+ Brown, G. E., Rev. mod. Phys., 31, 893. 
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and likewise for G-. The S’s are replaced by perturbation series like (3.7), 
with the result 


pes | dy... dt, T(b(e,t) pr(e', V(t) ««-V(tn))|0) 


n! 


G(reterst)=— 


(0) SY ay. dt T (VE) .. Vit,))| 0) 
<hr 


This differs from (3.8) in the appearance of y% and * in the numerator 
rather than V(0), and by time integrals from — oo to + oo rather than from 
—oto0. Evaluation by Wick’s theorem leads to diagrammatic represent- 
ation as before. The numerator diagrams now have external lines ter- 
minating at points (r, ¢) and (r’, ¢’); these represent factor pairings 
involving the extra /and:s+. Numerator diagrams with parts not connected 
to the external lines precisely cancel against the denominator, leaving a 
linked diagram expansion for G. 

The diagrams for G are exactly the same as those for Gin§3.7. The only 
difference is in the extended range of the associated time integrals. The 
quantity G, for negative t and t’, is actually given by 

wet gy SOL LM, t), or’, #))10) 
Geter oe <0 [Oy Ls, ste ek DO) 
The appearance of the model state (0| on the left replaces the quantity 
S(co, t) of (3.154) by S(0, t), and so restricts the time integrations as noted. 
We again obtain equations like (3.97) 
G=G-—iG?MG 
= G9 -$G MG?) iy ee te ae oe 


where M, like J, arises from proper self-energy insertions only. Unlike 
G and M, G and M are translationally invariant. They have therefore 
Fourier decompositions of the typet 


G(r,t;7r’,7') -|= ¥G,(E) exp [ik. (r—r’) —id(t—-t’)]. 


From (3.157) 
G, (6) =G,°(E 1 + M,(6) G,(E)}> 


2 G,(é)= {6- (or rie) a mi(s)} . . (3.158) 


Now in the first order of perturbation theory M;,, reduces to the usual 
Hartree-Fock potential, which depends on k but not &. Equation (3.158) 
then shows that the propagation is modified only in that the velocity 


} In order that Fourier transforms with respect to time exist in an unsophisti- 
cated sense, we again introduce some artificial damping. In perturbation 
theory this is conveniently done by adding (—ie) to the kinetic energy of states 
above the fermi surface, and +ie to that of states below. This corres- 
ponds to the modification (3-148) applied to the unperturbed propagator @®. 
The damping of G and M, constructed by perturbation series involving G° 
is then more complicated. : 
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dependent effective potential is added to the kinetic energy. More 
generally M can still be regarded as the effective single body potential 
but (because of the & dependence of its transform) it is an integral operator 
(non-local) with respect to time as well as space. In the theory of the 
optical model of finite nuclei (Bell and Squires 1959, Bell 1960) a quantity 
analogous to M is shown to give correctly the elastic scattering amplitude. 
The construction of this quantity is formally identical with that of M and 
can plausibly be approximated by the latter in the central region of a large 
nucleus. 

_ A natural generalization of the single particle energy of Hartree-Fock 
theory is given by the root of the equation 


/ ke? / 
Cr = 3M a M,. (€, e . ° . . . (3.159) 


This must be solved by iteration to yield an explicit series. The quantity 
é; is in general complex, indicating the transient nature of the particle or 
hole motion, which dissipates into more complex excitations. In the 
first order in v or in the K matrix, e;,’ of (3.159) is the same as the single 
particle energy of Brueckner theory §3.7._ But conceptually it is different, 
for it refers to real rather than virtual particles and holes; moreover it has 
meaning to all orders of perturbation theory, and not just in a particular 
approximation scheme. 

In the particular case k =k, the single particle energy is real and equal to 
the separation energy e, of §3.9. This is the theorem of Hugenholtz and 
Van Hove (1958). To see it, consider for a moment the case of finite ©. 
Then the propagator G,,(&) has a discrete set of poles whose positions give 
the energies of the (A +1) particle systems. Constructing from (3.159) 
the series for the energy level continuous with e,’ =k,?/2M in perturbation 
theory, setting «= 0 and then letting Q+0o, we must obtain a series for 
the quantity e,’. Moreover, this series is essentially the derivative with 
respect to k, of the energy per particle (§3.9). Now the series obtained from 
(3.159) for the optical model energy e,,’ in the limit k > k, can differ from the 
series just described only in the presence of the infinitesimal imaginary 
parts ze in energy denominators (for in this case the limit Q— oo is taken 
first). Thus e,’ as kk, is the derivature of a series which is obtained from 
that for H/A by inserting such imaginary quantities. However, because 
vanishing energy denominators do not occur there with finite weight, this 
makes no difference when the limit «>0is taken. Thus e;,’-e, ask—>k,. 

This provides, in conjunction with (3.137) another way of calculating the 
total energy per particle. It was so used for example by Galitskii (1958). 
There are other ways of obtaining the energy from the one-body propagator 
(3.146) or from the similarly defined two-body propagator, e.g. Klein and 
Prange (1958), Galitskii (1958). Such formalisms have the advantage, 
over that developed in § 3.3, of not introducing explicitly an origin of time 
and so working always with translationally invariant quantities. On the 
other hand, the total energy shift is not given directly but only after 
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for example an additional integration over ky, or over a coupling constant, 
or after a separate calculation of kinetic and potential energies. 


3.11. Final Remarks 

We have freely made rearrangements in the diagram expansion. In 
the next chapter it will be seen that summations over certain classes of 
diagram diverge for arbitrarily small interaction strength g for sufficiently 
large volume Q. The complete series may itself have this character, 
although that has not been proved. Despite this a rearranged expansion 
might remain valid. Suppose for example that (for fixed large Q) the 
rearrangement could be justified for small enough g, and that the ground- 
state energy varies analytically as g is increased to the required value. 
Then if the rearranged series also gives a result analytic in g it is the correct 
result. However, the ground-state energy may not vary smoothly with g ; 
for example, some originally higher level may come down and take over 
from what was originally the ground state. The considerations of §4 
suggest that this may in fact occur for nuclear matter. Although the 
finite temperature version of the formalism might then remain applicable 
(§ 4.4) the Brueckner—Goldstone approach of this chapter would no longer 
answer the right question. Indeed the particular rearrangements used 
may not provide a sensible answer to any question, for we will see that 
formidable difficulties present themselves in connection with them. 


§ 4. SINGULARITIES AND SUPERFLUIDITY 
4.1. Introduction 

The manipulations of the previous chapter have been made without 
regard to whether the series considered are suitably convergent, or whether 
the integral equations obtained actually have non-singular solutions. In 
fact, however, the Brueckner K matrix, expressed asaseriesin Vis divergent 
in an infinite system however weak the potential. Further, for nuclear 
matter there is good reason to believe that accurate evaluation of the 
results of the theory would lead to a singular expression for the energy. 
The finiteness of the results obtained in the previous chapter must then be 
attributed to some of the cruder approximations of the working (e.g. 
the coarse mesh used in the numerical integrations and the treatment of 
‘off energy shell’ effects in the Brueckner-Gammel calculations, and the 
use of only second order in the long-range part of the potential in the work 
of Moszkowski and Scott). These matters are discussed in § 4.2. 

In seeking some understanding of the singularity, we discuss in § 4.3 the 
effect of including additional diagrams in K, and in § 4.4 the generalization 
of the theory to non-zero temperatures. Finally, in the last four sections 
is considered an entirely different, essentially non-perturbative, approach 
to the theory, which leads to the suggestion that nuclear matter is superfluid. 
The basic theory for zero temperature is outlined in §4.5, and for finite 
temperature in $4.6. Some possible corrections are discussed in § 4.7, 
and § 4.8 contains explicit application to nuclear matter. 
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4.2. The K-matrix Singularity 


To illustrate this singularity (Cooper 1956 is the first published referencey) 
we use a separable potential, which we take to have the form 


{k’ |vo[k)=QO-19v,,0, id sade coro gpa (09) 


where k = }(k, — ky) and k’ = 4(k,’ —k,’), k, and k, being the initial momenta 
of the interacting particles and k,’ and k,’ their final momenta. We assume 
v; to be real and positive. Such a potential acts only ins states and so gives 
saturation if g is negative. The K matrix, in the effective field approxi- 
mation of § 3.7, for nucleons of total momentum p, is given from (3.17) by 


Ck’ | K( a k)=27gr,0, 
Q-1 One Vref (K", p)<k" | K(p)| k) 
g ” " 4.2 
2. e(k+3p) + e(k — 3p) —e(k” + $p) — e(k” — 3p) ae 
where f is defined in (3.60). The parameter & has been given the ‘on the 
energy shell’ value e(k+4p)+e(k—4p) since this is what is required in 


(3.122). We have used the notation e(k) for e,. Equation (4.2) is readily 
solved to give 


; 2-900 
(k’ | K(p)|k) = —— 20 Pe 4.3 
| K(p)|k) Tba5Uc p) (4.3) 
where 
1 Ope fk", p) 
be eee eI Pa 
2M P)= 6 2% ek” Ip) telk”— dp) —alk-+ Ip) -ak~ Ep) 
or, for Q+ o, 
dk” "2 
Stk.p)= | : oP eo SS Se 
Ik’ +4p1>ky (277)? e(k” + 3p) + e(k” — 3p) —e(k + 3p) — e(k — 3p) 
| (4.4) 


If e(l) <e(m) for all (>k,) and all m(<k,), then the integrand in (4.4) is 
positive. We can then obtain a lower bound for > by replacing the 
condition k,<|k” + $p| by the stronger one 

kp<k" +3p<ky+a 


where A is any positive quantity (which we later take less thank,). Fora 
reasonable form for e(k) we have 


Agrk” (k” — k) One A 
50 e(k” + dp) +e(k” — dp) —e(k + 4p) —e(k— 3p = 
(27)® e(k” + $p) + e(k” — dp) —e(k + 3p) — e(k — 3P) 
for all k” in the above range and for |k + 4p|<,, where A is some positive 
number. Then . 


td Be adie kp—-4p+aA—k 
;, =Al a EO RIES Coe fe ea (AG 
a AS aap ie kp+3p—k ey 


(4.5) 


whence >> 0 as k>k,, p>0. 
ie se ee ee 
+ Van Hove (1959) gives many references. 
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Asa first consequence of this note that, whatever the sign and magnitude 
of g, there is a region where > > 1; there the series expansion of K in powers 
of g is divergent. Moreover the associated power series for the energy 
shift is also divergent. This is of the form 


Oe = gS (—9)" An Po PEN sl Ae en ead 
where 
= dk__dp n 
Baa 12[ cay ay oP oe CP Rate fe) 


Imposing the more restrictive condition 


kp—-A<k+3p<ky 
and using (4.6), we find 
PA 302 8 


where B is a positive number independent of n, and 


Bret aaeet| cee eee A : 
=A 1 (yQ) —_________ - 
wtf def at | eee 


By change of variables and partial integration 
1 pl 1\ 
pas | die | dw (108) 
0 t %, 
1 1\” 
te | die Ce *) 
0 t 


=nP,_4—n!/4r+1 
whence 


n! 
P,=n!(Py- 1/12) +3 : 
so that, for large n, P,,/P,_,>n, and 


An|An-1 ZN. 


The series (4.7) thus diverges for all g. 

The divergence of the power series, even for arbitrarily weak non-singular 
potentials, is unexpected. But it is not a vital matter for Bruckner theory, 
whose final equations make no reference to such an expansion. A much 
more alarming consequence of (4.7) is that for g < 0 the denominator in (4.3) 
can vanish for values of k and p required in the evaluation of AH. Sucha 
zero occurs unless g is larger than the maximum of Y-!. The Brueckner 
expression for AH is then a meaningless singular integral. For a large 
but finite system the singularity may not occur for the allowed discrete 
values of k and p, depending critically on the value of g and the size Q of the 
system; as © org are varied the energy shift AZ oscillates infinitely. 
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The case of more general potentials} has been investigated by Emery 
(1959, 1960) again using the effective field approximation and considering 
only p=0. The ‘on energy shell’ K-matrix equation becomes 

K|m)=v|m)+y 21K m) keg liy 
yr 2e(m) — 2e(l) 
aN studied the homogeneous eigenvalue equation corresponding to 


(4.9) 


_— Hye] |m) 
NOs 52 e(m) —e(1) 


A necessary and sufficient condition for the series expansion of K, in powers 
of v, to exist is that none of the eigenvalues p are less than 1 in magnitude. 
Further, (4.9) has a non-singular solution if, and only if, there is no eigen- 
value of (4.10) equal to unity. Emery (1959) showed that as mk, there 
is always one eigenvalue that tends to zero, thus demonstrating that there is 
a range of m for which the series expansion of K is divergent. He aiso 
showed that with the central s-wave part of the Brueckner-Gammel 
potential (§1.3), and with the effective mass approximation (§3.7) for the 
e(k), there is not a non-singular solution of (4.9) in the triplet state if 
M*> 0-75 M nor in the singlet state if M*>0-:56 M. In the calculations of 
Brueckner and Gammel the effective mass at the fermi surface 


Wri. ch de(k) 


: Pere rm ae ans 


was found to be 0:73 M. 

Emery’s result for the singlet s-state was confirmed by Nuttall et al. t (1958, 
unpublished) who calculated the K matrix directly from the equations 
used by Brueckner and Gammel (§3.8). With A (see below (3.128)) chosen 
so that é(2e,,, 1) was continuous with e,, as m and I tend to k; (this makes the 
energy denominators come to zero, as in the effective field approximation) 
they found that a singularity occurred in (m | K | m) for m equal to 1-486 f~ 
(with k;=1-5f—). When A was increased the particle energies @(2e,, 1) 
just above the fermi surface increased so that there was an energy gap in the 
spectrum. The singularity disappeared when this gap became equal to 
about 0-4 Mev. 

Tn the calculations of Brueckner and Gammel the singularity would not 
have been noticed because of the coarseness of the mesh used in the 
numerical integrations. Also the value of A used in the calculations was 
probably such as to prevent the singularity occurring. The removal of the 
singularity by this crude treatment of ‘ off energy shell’ effects is of course 
unsatisfactory. Without this approximation there is continuity in e,, and 
el, 2e,,) at the fermi surface. 

The question arises as to whether these peculiar features of Brueckner 
theory imply some peculiarity of the system, or merely show that an 


+See also Van Hove (1959, 1960). 
{ Clement, C., Nuttall, J., and Sampanthar, S 
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unfortunate selection of terms has been made. For example, is the com- 
plete perturbation series (as distinct from the Brueckner series) always 
divergent? That this need not be so has been illustrated by Bell (1959) in 
the many-body problem with one-body forces. There one supposes the 
fermions to interact not with one another but with an external potential. 
A K-matrix approximation to the energy again leads to a singular integral, 
for a weak separable attraction, and a divergent series. However, this 
simple problem can easily be solved exactly (Fukuda and Newton 1956, 
de Witt 1956); the energy shift is an integral of an expression linear in the 
scattering phase shifts, having a convergent power series for sufficiently 
weak interaction. It has by no means been proved that this is not the case 
for the two-body force also. However, as we shall see in § 4.3, extra terms 
corresponding to those which remove the singularity in the one-body force 
problem do not do so in the two-body force problem. It is now generally 
believed that the divergence of the AK matrix is associated with the inade- 
quacy of the unperturbed fermi gas as a point of departure, and to the 
existence of a ‘superfluid’ state (§ 4.5). 


4.3. The Inclusion of Additional Terms in the Perturbation Series 


In § 3.4 we saw that it was possible to express the sum of a certain class of 
diagrams, the particle—particle interactions, as the solution of an integral 
equation, and in §3.7 we saw that certain additional diagrams could be 
included by modifying the energy denominators. We shall now show how 
to include still further diagrams which restore symmetry with respect to the 
fermi surface. Since the singularity discussed in § 4.2 is associated with the 
interaction of particles and holes near the fermi surface, this might be 
expected to provide some insight into the problem. 

The simplest extension of § 3.4 is to include all hole—hole interactions as 
well as all particle—particle interactions, i.e. to sum all time ordered 
diagrams of the types 


(4.11) 


Renney ric 4 2 ce ea aaens (EEL) 


with an arbitrary number of particle-particle and hole-hole interactions. 


e 
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This can be done (Chisholm and Squires 1958) by replacing (3.39) by the two 
equations 


{m,m, |K | m,m,)= (mm, |v | m,m, ) 


(mm, |v|h1,)<h1, | K | m,m,) 


+ ¢ 
on CRE IT ge 
CLE |K | m,m,)= «Il, |v] mm, ) 
> Hl Jo | hb’) h’1,’ |.« |m,m,) 
I, 4. (m7 + m,? —1,'2 —1,'*)/2M 
LE | | m,’m,’ )(m,'m,’ |v | m,m, ) 
+ Me Ee 4 
na, (m,'2 + mo? — 1,2 —1,*)/2M oo 


If we denote a K interaction by a heavy dotted line in a time ordered 
diagram, this last equation can be written diagramatically as 


one ee, elke baal’) 


+- 
+ 


which can obviously be iterated to give the desired result. For a separable 
attractive potential the additional terms included have the same sign as the 
particle—particle interactions of the Brueckner K matrix and consequently 
the power series is still divergent. With the same potential it was shown 
by Chisholm and Squires (1959) and Mehta (1959) that for a finite system 
the energy obtained from (4.12), (4.13) and (3.41)—with K replaced by K— 
again oscillates infinitely as the volume is varied (for sufficiently large 
volume) and has no limit as Q>oo. Thus, although an unambiguous 
solution, containing a 5-function, exists for the integral equation obtained 
by taking the infinite limit in (4.13), it is not possible to attach any meaning 
to the energy resulting from this ‘solution’. 

A further extension can be made by including in the sum. all ‘pair 
production’ diagrams, e.g. 


(4.15) 


: 


Indeed, if we use Feynmann diagrams (§3.3) these are not distinguished 
from the particle-particle and hole-hole interaction diagrams, since 
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the order of the interactions is not significant. Mehta (1960) calls all 
diagrams including hole-hole, particle—particle and pair production 
interactions ‘generalized ladder diagrams’. In §3.7 we saw that the 
contribution of all generalized ladder diagrams could be obtained by 


defining the operator K’ by 
K'=0-ikK' GGoy Wie a> eee 


where we are using the matrix notation of §3.7. For G@ we shall use the 
effective field approximation (3.119); then 


G,(t, t’)=exp[—ie(k)(t-t’)], © k>k, and t># 
= —exp[—ie(k)(t—t’)], k<k, and t<t’ 
=. otherwise,» <*. 6. Sie eae ee eee 


We can then define 
0 
K(6)= | exp (—tét) K'(0s4)dt. se) See 


—-o 


and, from (3.100), (4.16) and (4.17), obtain 


ap ake K'(é) Al) CLE |e 
le ee = e(1,) — (la) 8 
a y [Kms =Cm, es K'(é))| mym,){m,m,|v- (4 18) 
m,, mM, e(™my) + e(mg) — e 


The energy shift is given, from (3.105), by 


2 2 
AB=3 X (mym,[K’ (MEM) [mam (4.19) 

Equation (4.18) was obtained by Prange (1960). The contribution of 
the generalized ladder diagrams to the energy had previously been obtained, 
in a different form, by Galitskii (1958). A very simple method of obtaining 
the above result, due to Mehta (1960), is given below. 

Prange (1960) solved the integral eqn. (4.18) for a pure hard core potential 
in nuclear matter and found that the effect of the additional terms increased 
the energy shift by 10°% over that found by the Brueckner method. This 
is a surprisingly large correction to arise from the short-range hard-core 
potential, but confirms the result of § 3.5 that the low density expansion for 
the pure hard core converges slowly. No treatment of the above integral 
equation for K’ with a realistic potential is available, but the estimate of 
the hole-hole scattering correction in lowest order discussed in § 3.6, suggests 
that some cancellation occurs and that it is actually much smaller than 
with the hard core alone. 

However, for an attractive separable potential, the additional terms 
included in K’ have the same sign as those already included in Brueckner’s 
K, so the series for K’ is certainly divergent whenever the series for K is 
divergent. We can show further that, for such a potential, the energy 
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obtained from (4. 18) for a finite system becomes complex as the volume of 
the system becomes sufficiently large. It is clear, then, that the intro- 
_ duction of generalized ladders} does little to resolve the difficulties discussed 
in § 4.2. 

To obtain the above result, note that with the potential (4.1) the solution 
to (4.18) can be obtained in the form 


<k' | K(E)|k)=Q-g vy, vyF(@) «sw (4.20) 
where, for p=0, 


2B ( & 2 
F(é)= -$34 - ao Um ee )] (4.21) 
or 
2) — Igy uP BS Wd wm F (205) 
i ley jit oot O~ de, aa| ¢ OQ “de eee OE 
ye Sen (4299) 


Now this does not determine ¥(&), for clearly the quantities F (2e,,) can 
be specified arbitrarily. To complete the determination we specify that 
the solution be continuous with the iterative solution valid for small g 
(the manner of derivation obviously requires this solution). Now, from 
(4.21), one sees that the iterative solution has, for small g, singularities near 
& =2e, but not near &=2e,,. Suppose the quantities A; are those zeros of 


Z(é)=1+ 2 ae (4.23) 


as é > per, ier, 
each of which approaches an e,, as g—>0. Then, for F (A;) to be finite, we 
need, from (4.22), 
Um" F (2&m) 
SADR Ty See yea are Sea aed © SP 
2 a ereaea: 2€m — Aj ( ) 
for all i. This determines the A(2e,,). To find the energy shift we 
multiply (4.24) by the common denominator and regard the result as a 
polynomial equation for the A; The sum of its roots is then given by 


DA= Bn — GD en F (Zen). Oe See 4395) 


From (4.19), (4.20) and (4.25) the p=0 contribution to the energy shift is, 
ignoring trivial factors, 


AE, See RRL (ee )|m) 


Shy (Be, Agate ints = eee (4:26) 


m 


i ee ee ee ee eee 
+ The difficulties with the generalized ladders appear to have been observed 
first by J. Goldstone (unpublished). 
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This gives a reasonable result for sufficiently small |g] and Q. However, it 
is easily seen that on increasing |g| (with g<0) or Q, two of the roots of 
Z(&)=0, including the largest of the A,, become complex. It can easily be 
shown that the imaginary part of the largest \; becomes independent of Q 
as Q-> co, so that, since the other ), are all 0(1/Q), the imaginary part of the 
energy shift (4.26) is of the same order in Q as the real part. 

Mehta (1960) has given an instructive alternative derivation of the 
above results, exploiting the triviality of the one-body force problem. 
Consider such a problem with occupied levels of energy 2e,,, unoccupied 
levels with energy 2¢,, and an external perturbaing potential with matrix 
elements given by (4.1). The ground state energy shift is 


AB =D (Ag — 26h). spon ce one es eemeeeeaee 
where the 2’s are the lowest eigenvalues of the one-body Schrédinger 
equation 

S £2em Bact Kk o|k’y} Kk | D=ACk |p)... (4.28) 


These eigenvalues are in fact the roots of 


1+ Uae Ay Soe erg ee oa 
> ye 37 (4.29) 
Now, in terms of perturbation diagrams the complete solution of this 


problem is given by the analogues of the generalized ladders (there are no 
others). A typical such diagram is 


Pre a ev Ct), 


The sum of all such diagrams must give (4.27) and (4.29) exactly, as can be 
verified. The two-body force diagram analogous to (4.30) is 


oe 
ae 6481 
. 
| 


To obtain the contribution arising from (4.31) from that arising from (4.30) 
we recall the rules of $3.3 and introduce ‘ ya Hoe (—1) to account for the 
extra closed loop and (ii) a factor (— 1) for each backward line in (4.30), to 
account for the fact that they appear twice in (4.31). The latter can be 
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effected simply by replacing v,,2 by —v,,?2 in all expressions. Thus, from 
(4.27) and (4.29) we immediately obtain the result (4.23) and (4.26). It 
will be noted that the difficulty of complex eigenvalues never occurs in 
(4.29), and this is an essential difference between the one and two-body force 
problems. 

Finally in this section we recall that, in the derivation of the Brueckner 
approximation starting from a wavefunction with two-body correlations 
(§ 2.5) we obtained an expression which included the energy shift (m | K |m) 
in the energy denominators (see (2.43)). This could be obtained from the 
low density approximation of §3.4 by modifying the hole propagators to 
include the interaction of the hole with itself, i.e. 


; kia ooae ey) 


Gottfried (1958) considered the general solution of this modified K matrix 
equation for the separable potential. For simplicity we take total 
momentum zero and put (m|K|m)=A,, whence 


eae ee ee 
ee a” Opel 
Be ee Pie i intent Sayin Poi ee NS rg 5 
Q vm Q > ei — Em is > (e bie Cm) (ey a Con ae Ay) ‘ ( 


If the resulting equation has a solution with A, > 0 as Q— o it is reasonable 
to neglect the A,, in the denominator of the last term of (4.33). The 
quadratic equation for A,, then has solutions 


2 1/2 
Di 33 | Bn? ik 7 Am on | 


1 Pav) eT Lee es a Oe Ao ae( £54) 
ie 2gA m 
with 
ida Sse eee y (4.35) 
Ke QO I (eG)? 
eae ae FMS BN ait ON tae 
Boal + Ghee (4.36) 


Away from the neighbourhood of the zero in B the second term under the 
square root of (4.34) is small compared to the first for large Q and we can 
expand, thus obtaining the usual solution if we take the negative root, i.e. 


ee mpgs 8 os oe Sere eae as (A 87) 
m 


P.M.S. Yi 
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However, near B,, = 0, the solution of (4.34) does not have a pole but is given 
by 
v 

Palas Pod Qie 4 iB ss tUge ne CONES Reena ae 
This is 0(Q—/?) for large Q so the condition mentioned below (4.33) is 
satisfied by this solution. The usual form (4.37) is 0(Q~). The region 
over which the true solution differs appreciably from (4.37) is 0(Q7"?) so 
that ¥,,A,, is 0(1). 

This procedure therefore gives an unambiguous, finite answer. From 
the viewpoint of perturbation theory the inclusion of the modification 
(4.31) to the hole propagator, without, for example, a similar modification 
to the particle propagator, seems rather artificial, so it is difficult to ascribe 
much significance to this result. Note however that for a potential which 
acts only in states of total momentum zero (§ 4.5), only the hole propagator 
is modified in this way. For such a potential, with constant matrix 
elements, Katz, de Shalit and Talmi (1960) have shown that this treatment 
gives the exact answer in the strong coupling limit. However we believe 
this result to be largely fortuitous, since it is not true for example when the 
matrix elements of v are not constant. To see this we take the strong 
coupling limit (g—> oo) of (4.33) which gives 


An ait ( Cant + ~ vp ) . eS i aa ae (4.39) 
1 


We consider a situation (see § 4.5) where the number of momentum states 
k is P, each of which can be occupied by a particle with one or two spin 
values. If there are N particles and if the interaction is just between 
particles of total momentum zero and opposite spin, then 


AH=>A,=9 ( Sv,2+ L307 ) et ie re) 
m m 1 


This reduces to the exact answer (4.66) when v, is independent of k, but 
is quite implausible in general because of the discontinuity in the dependence 
on v,,and 7. Nevertheless the discussion has value in indicating with what 
delicacy terms of apparently negligible order in Q-! must be considered, and 
that the self consistent theoryt may then give results which are at least 
finite. 


4.4. Ladder Diagrams at Finite Temperatures 
It might be suspected that the difficulties discussed above are due to 
concentrating attention on one of a nearly degenerate set of unperturbed 
states, i.e. on the unperturbed ground state which is only separated from 
many excited states by an energy of 0(1/Q). It is possible that perturba- 
tion theory should work better when applied to the mean energy of a band 


of such states. To introduce averaging over many low excited states the 
oo PE eee ae SS antes Feats Ate So WS. Me EN 


} See Kimmel (1961), Wild (1960), Katz (1960 a). 


Theory of Nuclear Matter 291 


most obvious method is to consider the system at a small but non-zero 
temperature, and allow the system to become infinite before passing to the 
zero temperature limit. It turns out, however, that the integral equation 
giving the formal sum of ladder diagrams has no non-singular solution 
when the temperature is reduced below a certain critical value, independent 
of Q for large Q, which is, moreover, the same as the critical temperature 
obtained in the theory of superconductivity to be discussed in the next 
section. 

We first outline briefly the general formalism for non-zero temperature, 
following Matsubara (1955) and Bloch eau de Dominicis (1958). The 
grand partition function is defined by 


Z=Traceexp[—B(H—-AN)] . . .°. . (4.41) 


where H is the hamiltonian operator, V is here the number operator, A is 
the chemical potential and B = (k7')-!, k being Boltzmann’s constant and 7’ 
the temperature. We write tae in the form 


Z=Z, <S(B)) Bee ep keene 22 wetas. (4,42) 
where we use the notation 
seal ead mt eee 


and where Zp, py are respectively the grand partition function and the 
statistical operator in the unperturbed system, 


Zig WE GXD le Olt1 SAIN) Ye fais et ee ov (4.43) 
and 
pPo=Zy texp [—B(Hy)—AN)]. Seal een fur se (444) 


The quantity S(f) is given by 


S( Ge teXp (Gig \exDp (BH). tere. ae yn (445) 
and satisfies 
=~ ASE) concent = (4546) 
where 
V(g)=exp(BH,) Vexp(—8H), .... . «2 (4.47) 


V being the perturbation. 

Equation (4.46) is formally similar to the equation for the rate of change 
of the S matrix of field theory with time, and the methods used to simplify 
the perturbation expansion of the S matrix can be used here. The boundary 
condition is that Z=Z, at B=0. Thus 
ibe Seip 
Re aes 
where 7’ is the ‘time ordering’ operator which in this case rearranges the 
operators V so that they occur with decreasing values of uw when read from 
left to right. We now use a generalization of Wick’s theorem (see 


| CH. di CIV Gn) AY Manne ea 
0 


wes 
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Matsubara 1955, Bloch and de Dominicis 1958 and Gaudin 1960a), 
which enables us to express (4.48) as the sum of contributions of all 
Feynmann diagrams (§3.3) with no external lines. We can explicitly 
eliminate unlinked diagrams (cf. § 3.3, see Bloch and de Dominics 1955 and 
Thouless 1959 b for details) and the result is conveniently expressed in terms 
of the thermodynamic potential ./ defined by 


Z=exp(—f¥). re a set es et gee) 
We obtain 
1 2 (-1)"/(° 
AGB=A,--> I. du,...du,(T[V(u,)... V(u,)]>;, (4-50) 


n! 


n=1 
where, as before, the subscript L means that only connected diagrams are 
to be included, and where .Y, is the unperturbed thermodynamic potential. 
The contribution of a given diagram to (4.50) is obtained from the following 
factors : 

(1) A factor —4v(x, —y,,) for an interaction line from x,, to y,. 

(2) A factor G°(B; x, w’) for a line from an upper end 2’(=x’, wu’) to an 
upper end x, and similar factors @°(8; y, y’) ete. for a line from a 
lower end yz’ to a lower end y, etc. 

(3) A factor (— 1) for each closed loop. 

These factors must be integrated over all space and over values of uw from 
0 tof, and multiplied by (6n!)~'. The propagator G° is given, in momentum 
space, by 

_ exp ((2/2M)—A\u-u') 


G,(B;u, u') = —T+expA(@pMoay’ U >U, 
_ exp ((/2M)—Alu-w) oy (4.51) 


1+exp—B(k?/2M—)) ’ 
In the manner of § 3.7 it is possible to introduce an effective field modifica- 


tion to the propagators. We will not go into details but will allow for this 
possibility by replacing k?/2M by e;,. We shall use the notation 


e(k) =€,=e,—X. . 4 « = 5 : (4.52) 


Following Thouless (1960 b) we consider the sum of all generalized 
ladders, i.e. diagrams in which two lines enter and pass through all vertices 
in the diagram. A typical diagram of this kind is shown below (3.97). 
We denote by L,,,(K; wu’, uw) the sum of all such diagrams where the 
incoming lines have momenta K+m and the emerging lines momenta 
K+n. Then, considering the case where the particles have opposite spin, 
we have 


Diy pl K 5 0, U) = Gi 4 (85 U', UG _ (8; u’, u) 8, 2 
B 
> i dU es (BU, WV (8 5u',0")<n]v[k) 
0 
x Dim (K 5 u”, w) +) OT eT Rag wD, al ins RN Uhh) 
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It follows (Thouless 1960) from this equation that L depends on uw, and w’ 
only through (w—vw’) and that it is periodic in u—w’ with period 8. Thus 
we can write 


(eo) 


DB, A (K3 usu’) = NS a bps m(K, v)exp [2rvi(w — wu’) /B] (4.54) 
where rele io 
L,,..(K,v) = (amb (eBe(k + n)] + tanh [$8e(k—n)] 


; 2Ble(k+n)+ e(k—n)]—47iv 
x | bam —BE CO] Li (Ks 9) |. er een IC easy 
k 


This equation has a non-singular sclution only above a certain critical 
temperature. Assuming the singularity to appear first for k= v= 0, this 
temperature is the lowest temperature above which the homogeneous 
equation 
1 

va EUs, Ul atic ye ke aS) 
; > €k k 
has no non-trivial solution. This critical temperature is identical with 
that of the Bardeen—Cooper-—Schrieffer superconductivity theory (§ 4.6). 

Thouless shows that there is always a non-zero critical temperature for an 
attractive separable potential. More general potentials are considered by 
Emery (1960). He finds in particular that there is a non-zero critical 
temperature when the Brueckner K matrix is singular; the converse of this 
does not hold. 

The statistical perturbation theory reduces in a trivial way to the previous 
ground state formalism when the temperature is made zero for fixed Q. 
To see this it is essential to note that with finite Q at zero temperature the 
particle number is a discontinuous function of A and the interaction strength. 
The more interesting case, that the limit Q— oo is taken before 8B > 0, 
requires more discussion (Kohn and Luttinger 1960, Luttinger and Ward 
1960, Bloch 1960b). 


4.5. The Bardeen—Cooper—Schrieffer Theory 


The divergence of the K matrix arises because, for the interaction of 
particles with momenta sufficiently near to the fermi surface, the allowed 
energy denominators can be indefinitely small compared to the strength 
of the interaction. Such a situation corresponds to a ‘strong coupling’ 
problem, and suggests that further understanding might be obtained by 
consideration of the strong coupling limit. To this end we first consider a 
system with 2P degenerate single particle states |k); k= +k,, +kg,... + Kp, 
where the state +k has opposite spin and momentum to —k. We ignore 
isobaric spin in this example. Let the interaction be between pairs of 
particles with total momentum and total spin zero, and be otherwise 
independent of k, i.e. 


V=g 2, Oe Seas ea Maen) Feat = (4008) 
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where the sums are over k and k’ equal to k,, ky... kp. We put 


b= Oty bhatt te ee ae ee en 
then 


Vagy bythe, 2 2 (4.60) 


It is easy to see that the pairs of particles corresponding to the 6 operators 
obey Bose statistics when the momenta are different, i.e. 


(6, Opt} = [6,0] = [0 Oe =0sa tor (ks kee os) 
For k=k’ they obey fermi statistics provided the space of states is limited 
to that for which only pair occupation is allowed (i.e. if the state +k is 


occupied then — k is also occupied and vice versa). Thus, in this restricted 


space, 
[b, both eee ek eer 


Since the hamiltonian can be expressed entirely in terms of the b operators, 
the energy eigenstates are in the restricted space in which only pair occupa- 
tion is allowed. Asa simple example suppose there are only two particles 
in the system, then the basis states are b,*+|vac.) and, if we take their 
unperturbed energy (equal for all k by hypothesis) to be zero, we have for 
the matrix elements of the hamiltonian 


(KTH [kya g Sr ee ae 


This can easily be diagonalized; for an attractive potential all the eigen- 
values are zero except the lowest which is —Pg. The ground state is 
given by 


( be") [vos 
k 
By analogy we expect the ground state for NV particles to be given by 
N/2 
[0) = ( by) | vac. ). . . . . . . (4.64) 
k 


Operating on this with the hamiltonian we obtain 


N N 
1 [0)=9 (PF +1) 0); tiny See el See 
so that it is indeed an eigenstate, with energy 
_N N ; 
Ab=g(P-5 +1). ° . : . . (4.66) 


For a proof that it is the lowest eigenstate for g <0, and for a derivation of 
the complete spectrum we refer to the original papers (Wada et al. 1958, 
Mottleson 1958, Thouless 1960 a, b and Katz et al. 1960). 
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We now turn to the general problem where the single particle states are 
not degenerate, but we still restrict the interaction to be between pairs with 
equal and opposite momenta. With nucleons we can consider parallel or 
antiparallel pairings of isobaric as well as ordinary spin. For simplicity 
we take both antiparallel. Indicating ordinary spin explicitly and 
absorbing isospin into the symbol k we have then the ‘ reduced hamiltonian ’ 


eas 2 &( Ut Ayer + Aye y Tey) + >» (KOK) ar F aig Fay dyy re 

(4.67) 
The e;, are single particle energies which may contain effective field modifica - 
tions from parts of the interaction not now considered explicitly. 

The interaction retained in (4.67) is of order Q- relative to the total 
interaction. It is found however to cause an energy shift per particle 
independent of Q. It is hoped that certain qualitative features of the 
eigenstates of H®°S may remain true for the total hamiltonian. One 
cannot be entirely confident of this because the neglect of interaction in 
states with total momentum not zero gives the hamiltonian a peculiar 
character—in configuration space it has an infinite range of non-locality. 
Incidentally it should be noted that although H®°S allows the occupant of a 
given state to interact with that of only one other state, in perturbation 
theory it generates other diagrams in addition to the ladders. 

The natural generalization of (4.64) to take account of the different 
energies of the unperturbed states would be 


Np 
& Wy My +41?) | vacuum ) 
k 


with w, to be determined variationally. However it is more convenient to 
work with the state (Bardeen et al. 1957, referred to as BCS) : 


posy qt okt Cutt] 
[PBce = Cs ~ (+l PE? 
where the product is over all momenta and isobaric spins. The variation 
parameters g are related to the occupation probabilities h for individual 
SES Py = |e PCE + | I |): Gpneamaps> 2269) 
The state (4.68) is not an eigenstate of particle number, but turns out to be 
sharply peaked about the mean value; the projection on to a state of 
definite particle number has the required form. 

We must minimize the energy expectation value for given mean number 
of particles. Introducing the Lagrange multiplier A we find 


ree] {Bcs NI eee 


= 524 1a ae elaine 
k 


Silat sae ie aac 5 (0 
e+ lo on 
(4.71) 


where as before qe f mle 


Jvacuum) . . . (4.68) 


Jue 
sea )(1+] 9 |?) 
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The last term of (4.70) is of order Q- relative to the others and will be 
ignored. Variation of the g, yields the stationary conditions 


be Ne + 2€% Ik — Mae =O (4.72) 


where 
— klo[k’) ee ee re wes 
¥ dklo|k’) 
and where we have assumed (k|v|k’) to be hermitian. The solution of 
(4.72) is 


1 
t= [eet [oe ee) ETA) 


where we have taken the positive square root since, as we see by forming the 
second derivative of (4.70), the other root gives instability. Then, from 


(4.69), 
h =1[1- aoe . 4: ey op 
fee [ex | ex Py? ( 


This is approximately unity for «<< —j, and zero for e, > 4,, and varies 
smoothly between these limits, so that the sharp cut-off at the fermi surface 
for the unperturbed gas is smoothed out by the interaction (provided 4,40 
of course). 

If we put (4.74) into (4.73) we obtain 


(kok Dee 
=— ) a: if Se aig ce CEE 
lid © 2fe,? + [ue PPB ( ) 


One solution of this is the trivial one 4, =0, for which we obtain the energy 
and occupation probability of the unperturbed state. In certain cireum- 
stances this equation also has a non-trivial solution. If this solution leads 
to a lower ground-state energy than the unperturbed energy, it follows that 
the true ground-state energy, being at least as low as this, is lower than the 
unperturbed energy and therefore unattainable by perturbation theory 
(which gives a zero energy shift for H®°8), 

As a simple example which shows the above phenomena, BCS take a spin 
independent potential of the form 


(k’ |v [k)= = "0 for k,—-t<k<k,+t, 
ae ky—t<k'<ky+t, . .. . (4.77) 
= 0 otherwise, 


where t<k,. Then, from (4.73), ;4,=€y, say, is independent of k and 
(4.76) becomes 


Leu t 
POse 0 «oD (e,2-+ <>) +e) (4.78) 


where >’ means the sum over all momentum states with ky—-t<k< ky+t. 
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Obviously this has a non-trivial solution only if the potential is attractive, 
Vp > 9, when, in the limit of infinite volume, 


de 


sinh (1/n v9) en 


Eg = 


where Qn, is the density of momentum states per unit energy interval near 
the fermi surface, given by 
134A 
ea Selihe vit hates 24 t (4,80) 
k 
ke LF 


Qn,= Ee 
de=|[e(k,+t)—e(k,)| = |e(k,)—e(kp—t)]|. . . . (4.81) 


0k 


and 


To find the energy it is necessary to calculate A. Alternatively we can put 
A=e, (its value in the unperturbed state) and then calculate the number of 
particles being considered. With A=e, there follows from (4.70), (4.75) and 
_ (4.79) after some calculation 


(BCS | HBOS _)N | \PBOS)_ y'(e,, —)) 


# 2 
- 2n;(de)" (4.82) 
exp [2/n,v)]—1 
We first use this result to show that the system must have an exact 
ground-state energy which is lower then the unperturbed value, and hence 
unattainable by perturbation theory. To see this we note that, for any 


normalised Y, i" 
(¥ | HB°8_AN |'¥) > E(N)—-AN Care. (4.83) 


where L(V) is the exact ground state energy of H®° for NV particles, and V 

is the solution of 

ok(N) 
oN 

If the unperturbed energy is equal to the exact energy then the solution of 

(4.84) is the number of particles in the unperturbed state, i.e. 


Nias liee Ng ae cea oo hat kx,’ (4.86) 


(4.84) 


A=e,;= 


In this case (4.82) contradicts (4.83) with ‘V8°° for ‘Y. Since the unper- 
turbed energy cannot be lower than the exact energy the result follows. 

In general (4.85) is not satisfied so that (4.82) does not immediately give 
the variational approximation to the energy shift. It does, however, give 
it correctly to first order in the energy shift, since, to this order, 


E(N)~ E(N,) +N —N,). 
Then, from (4.83), 


CPBCS | HBCS_)N | PBCS) _ ¥(e,,—A) > H — 2m =AH(N 
(4.86) 
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When the BCS trial wavefunction is used for the strong coupling problem 
discussed at the beginning of this section (Katz et al. 1960), 4, and g, given 
by (4.74) and (4.76) become independent of k, and are thus determined by 
the condition that 


2 
hi; oops le Ee one ei aes ees 
2h 1+] 9? 
Substituting for g, in (4.68) yields on evaluation of the energy 
N NE hay, 
=g—(|P-—+—,; oA, «cual ile” oe eR 
AH =g 3 (P a =p) (4.88) 


which is to be compared with the exact answer in (4.66). The difference 
is due to the fact that the BCS trial state gives an average over a small 
region of particle numbers. The error becomes negligible as P— co 
(assuming P— N/2 is not 0(1)). 

In § 4.8 we shall investigate whether (4.76) has a non-trivial solution for 
nuclear forces. Here we shall assume such a solution exists and consider 
the possible properties of the system. To do this it is convenient to intro- 
duce the idea of ‘quasi-particles’ by means of the Bogoliubov—Valatin 
transformation (Bogoliubov 1947, 1958a, b and Valatin 1958). 

We put 

ae As = Jun y* 
ste di 


é — Gey THe Ue * 
ky “T1+] 9 ad eo ; 


It is easily seen that for any g, the operators €& and &, obey the fermion 
commutation relations and that (4.89) is a canonical transformation. The 
inverse transformation is given by 


(4.89) 


_ Sur tone Sry 
eT T+ [ge PP 
iy Sk Gee ke 

[1+ | gx [PJM 


If the g, are identified with those in (4.68) the state |‘YB°S) is the ‘ quasi- 
particle vacuum’, i.e. 


(4.90) 


AK 


Eg | PECs 0. hoe tea) ee aa 


In the BCS theory the quasi-particles are regarded as independent 
particles, in which case the excited states of the system are given by 


Sige Ciao? ee Cine t Hegre ; n=l, Qyeve ee ee (4.92) 


Various attempts (see below) have been made to justify the use of these 
states as the excited states—they are not of course exact eigenstates of 
H®°S: the real justification, however, lies in the success the theory has had 
in providing a basis for the theory of superconductors. To find the ‘ energy’ 
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of the states (4.92), ie. the expectation value of the hamiltonian in the 
States, we insert the transformation (4.89) into (4.67). Then, ignoring the 
interaction between quasi-particle, which we can do if n is small, i.e. 0(1) 
rather than 0(Q) for large Q, the part of the hamiltonian which is diagonal 
in the space of states (4.92) is 


bol 


9 (k |v | k’) 
(HBOS _ : pea get: os Due Ju « 
see 2 Te Tg Tt SO [oe loa 


€% Ce eae) (9u* Dre (Due Jue + Je Jue™) {Kk |v | k’) | 
+ a 
rear 1+ |g | 2 T+ [oe FI oe [1+ | ue PILL + [Qe [7 


tee y Me EIT GL, deme nl e e. o (408) 


Using (4.73) and (4.74) we find for the energy of the states (4.92) 


E({k3) = > Len? + [o,,[? 1? + BBS Net Pour ed G4) 


where H®°S is the ground-state energy in the BCS approximation. 

The lowest excitation energy will be obtained for a state with one quasi- 
particle with (for a reasonable potential) |k|=,. Its excitation energy 
will be p, z We can expand p, in partial waves (Brueckner, Soda, Anderson 


and Morel 1960, Anderson and Morel 1960) 
Pee a LAK Ka ap Fae ota nee (4. 9D) 


l,m 
where k, is some fixed vector, and suppose (see § 4.7) that the non-trivial 
solution occurs for just one partial wave. The particular ,/™ will then 
depend just on the interaction in the /th partial wave state. If this is an 
s-state it will lead to an energy gap in the spectrum of the complete system 
since the lowest energy will then be that of a particle-hole pair with energy 


2ey= 2H, . Ser an rey Ew ert 49G) 


For other values of / the energy spectrum of single quasi-particles will be 
given, in this approximation, by 


Dy = pig? XY” (Ky Kole ae es 1 (4.97) 


For suitable directions this will go to zero so there will not be an energy gap 
but only a reduction in the density of low lying states. For these cases the 
material will become anisotropic. Some discussion is given by Anderson 
and Morel (1960). 

' There are alternative ways of obtaining the results of this section. 
Bogoluiboy (1958 a) used the transformation (4.89) and chose the g, so that 
no zero energy denominators arose in a perturbation expansion starting 
from |¥®°S) as the unperturbed state. Bayman (1959) used a trial 
wavefunction with a fixed number of particles and evaluated the expectation 
value of H®°S by a saddle point approximation. 
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4.6. BCS Theory at Finite Temperatures 
We briefly outline the generalization of the theory to finite temperature, 
following BCS and Valatin (1958). It can be shown (Peierls 1938, Valatin 
1958) that it is necessary in this case to minimize the thermodynamic 


potential 
A = (HB°S—AN)—-TS twee ar a eet es ee tra 


where the bar denotes an ensemble average and S is the entropy. We 
calculate these quantities by assuming we have a fermi gas of quasi- 
particles, defined by (4.89). Then 


TS=—B dh log f+ (1 —fy) log (1-f.)} ee et, (4.99) 
where 
Fe= Err? bet = Seat Ene): Mie Se ee (4.100) 
Using (4.67) and (4.89) we obtain 
= ke Berta 1—2f,) 
+3 dklolk’) He (1 2f,)(1- fy) - TS. 


ea 1+]9eP 1+ [oe P 
(4.101) 


On minimizing with respect to g, as before we obtain (4.74) but with (4.76) 
replaced by 
(k [e[ k’ (1. = 2h) 


Se et eee ae 4.102 
Me a Shes 24 os jaye MK’ ( ) 
Minimization with respect to f, yields 
Sf, =f +expB {e,2+]u, PHAR, . 2 2... (4.103) 
whence (4.102) becomes 
k |v |k’ - 
[y=-> See tanh [$8 (€,7 + |py-|?)"?] py (4.104) 


wr 2[eg* + | ee PR)? 
This can be written in a more convenient form by introducing the definition 
_ tanh [$8 (€,?+ | ex [?)!?] 


a] Ze. + | pn [2 a Sa 
giving 
tanh [48 ( Oana ih oh : 
X= - i ic nel PE oo (k vlk Sy Pe 4.106 
lett | pr a7 > | I Xe ( ) 
The critical temperature 7’, = 1/k8, is the highest temperature below which: 


this equation has a non- trivial solution; above this temperature the 
thermodynamics is that of an ideal fermi gas. We can obtain 7’, more 
simply from the linear eigenvalue equation 


tanh ($fe,,) 


Aa on DTep! 


> (klolk’> x. oro tae kre 
= 
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By arguments like those of § 4.8 the existence of an eigenvalue in the range 
0<p <1 is equivalent to the existence of a non-trivial solution of (4.106). 
For all eigenvalues (Emery 1960) || increases monotonically to infinity 
with increase of temperature. Thus 7’, is the highest temperature at 
which (4.107) has an eigenvalue .=1. Below this temperature there is 
always a non-trivial solution of (4.106), and nowhere above it. Comparing 
(4.107) with (4.56), 7’, is also the temperature below which the sum of 
generalized ladders is divergent. It follows (Emery 1960) that when the 
Brueckner K matrix is singular the BCS theory gives a supernormal state 
at zero temperature. 

Bogoliubov et al. (1958) have obtained these results by dividing H®°* 
into two parts, one of which is exactly soluble. For two different choices 
of the soluble part the effect of the remainder is of relative order Q-! in 
every order of perturbation theory. These give respectively the unper- 
turbed and BCS free energies, and it was asserted that both solutions are 
exact for a large system. But such reasoning is clearly not conclusive (see 
Wentzel 1960). On rigorous examination (Bogoliubov 1960) the unper- 
turbed solution has been rejected, and the asymptotic exactness of the 
BCS solution reaffirmed. 

Gaudin (1960 b) has found a part of the diagram expansion in statistical 
perturbation theory which yields the BCS result on formal summation. 
It has also been obtained by Martin and Kadanoff (1961), following 
Gorkov (1958), with the Green function method. 


4.7. The Effect of the Residual Terms in the Hamiltonian 


Here we consider briefly the effects of H—-H®°%. Although this is the 
main part of the potential it is hoped—and this is one of the principal 
assumptions of the BCS theory of superconductivity—that its inclusion will 
not qualitatively upset the main features of the above model. 

First of all it is necessary to decide which particular spin and isobaric 
spin state should be considered in the interaction retained in HP°. It 
seems reasonable to include that corresponding to the ‘strongest’ attrac- 
tion, where this might be defined as that which gives the largest critical 
temperature. 

Thouless (1960b) makes the quasi-particle transformation (4.89) and 
then considers H — H ®°S as a perturbation on the quasi-particle vacuum. 
He finds that the energy gap and smoothing out of the fermi surface 
associated with the previous solutions of H ®°S are sufficient to make the 
sum of quasi-particle ladder diagrams, with total momentum different 
from zero and spin equal to that of the interaction retained in H?®, 
finite at all temperatures. In fact this sum goes as logp, where /p is the 
total momentum, for small »(40). For p=0 the sum is divergent as 
previously stated. The sum of ladder diagram with different total spins 
is finite at the critical temperature but not necessarily at all lower tempera- 
tures. Thouless (1960 b) suggests that if it becomes divergent it probably 
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a 


implies the formation of ‘bound clusters’ of more than two particles. We 
know of no investigation of this problem appropriate to nuclear matter. 
Anderson (1958) has shown that the longitudinal oscillations, generated 


by the operator 
Pa on Deer A, OT en A bd (4.109) 
k 


acting on |‘ ®°S), have energy expectation values (of H) continuously 
distributed above the BCS ground state energy. Thus the energy gap 
is removed by these collective excitations, although the density of states 
remains much reduced. For metals the coulomb interaction has the 
effect of raising the energy of these collective state so that they do not fill 
the gap. The collective states are important in considerations of the 
electromagnetic properties of superconductors. For these and other 
details relevant to superconductors we refer to the original papers mentioned 
above and also to Kuper (1959), Abrikosov and Khalatnikovy (1959), 
Bogoliubov et al. (1959), Rickayzen (1959), Nambu (1960) and Wentzel 
(1960) where further references are given. 

As we have stated the main justification for the BCS model lies in its 
success in treating superconductivity in metals and therefore, since the 
physical situation in nuclei is not the same as that in metals, we may not be 
justified in applying the model to nuclear matter. In particular we have 
to remember that the hard core demands considerable further correlation 
beyond that contained in ‘VY ®°S, Further it is known that self-consistency 
effects are important in nuclei and the BCS theory gives no prescription for 
dealing with such effects. In view of these points, and the fact (see §§ 4.2, 
4.8) that the existence or not of a non-trivial solution of the BCS integral 
equation depends sensitively upon the treatment of self-consistency, it 
seems premature to make confident assertions that nuclear matter is, or is 
not, a superfluid. 


4.8. The Superfluidity of Nuclear Matter 
With the above qualifications, we now discuss whether the forces in 
nuclear matter are such as to cause superfluidity{. To this end we could 
consider solutions of (4.76); alternatively we can argue as follows. We 
introduce 


Yk 
x eae ais een ae es 


Erato eat aloe 4.110 
V+ |g l? — 2fe,? +] oe PIE? ! 


then, using (4.70), 
(pBes [HBCS Sry, \¢ | wBCS) _ > (Cm —A)= >| ee \{ Fadtie 4 | Xk [2)12} 
k 


+4 D Ck [op < yxy (4-111) 
k, k’ 
Then superfluidity occurs whenever there exists a x which makes (4.111) 


+ Aprikosovy, A. A., and Kuauatnikoy, I. M.. Advance. Phys., 8, 45. 
{See also de Dominicis and Martin 1958. 
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less than zero. A necessary and sufficient condition for this is that there 
should exist a y, such that 


ZL xP lex l+ 3 & xe xe (ke |v | k’) <0. Air e(4 sto) 


This is necessary because the latter expression is less than the former, and 
sufficient because they differ negligibly for small y. Finally, the above 
criterion is equivalent to the existence of a negative eigenvalue é of the 
equation 


Jel xe + BKK Xe = Exe Pa co (4-112) 


a condition which was obtained by Bogoliubov (1958 b) directly from the 
condition that the normal solution should be unstable. 

This last equation was used by Solovier (1958) who allowed for the 
possibility of different spin and isobaric spin states being responsible for 
superfluidity. His equations are only valid for weak potentials, and he 
appears to have required that there should be a superfluid solution in all 
angular momentum states in which the interaction is important. Thus he 
obtained several conditions which have to be satisfied by the potential. 
In fact, for a sufficiently weak potential, it is necessary and sufficient that 
{k;[v|[k,) should be negative for one particular spin, isobaric spin and 
angular momentum state in order that (4.76) should have a non-trivial 
solution (tensor forces prevent these being good eigenstates in general and 
obviously introduce further complication). 

Cooper et al. (1959) used (4.112) as the basis of a variational calculation. 
They considered interaction in s states for which, with 


MiMi) eV met cer yore Sg wt ou (Seb LD) 
the appropriate condition is 
2[ “Lexis Pae+ | “or enfear <0 va (48118) 
0 0 
where 
fle) =rxle)=r(2n)-% | exp (ik F) xy Hen (Aelia) 
or, from (4.115), 
2 1/2 fo : 
{n= (=) [fee sin ery. Pe de een 4 ELS) 
7 0 


For e, they used the effective mass approximation; presumably with the 
value found by Brueckner and Gammel for the effective mass at the fermi 
surface. With the trial functions they used they were not able to satisfy 
(4.116). Mills e¢ al. (1959) continued this work with improved trial func- 
tions. For e;, they considered three possibilities. 
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(i) An approximate polynomial fit to the Brueckner—Gammel (1958) 
values (presumably the curve e, in our fig. 2, though this is not made clear) : 


€, = ?/2M +[—111+ 39-1 (k/k,)? — 8-26(k/k,)?] Mev, for k<2-4k, 
=0 for) k< 24k; 


(ii) Re Vo 

ey = om sf ieee set ah Paes (4.119) 
with V, and « fitted by the requirements that e;,, should be the separation 
energy and the effective mass at the fermi surface should be the Brueckner— 
Gammel value (which they incorrectly take to be 0-671), and 


(iii) €,=k?/2M*, M*=0-67 M, 
i.e. the effective mass approximation. 


Fig. 2 


ry i 
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50 
0 
Beis ke) 
-100 
0 0.4 0-8 l-2 1.6 2-0 2-4 
k/k 


Typical form of the self consistent energy spectrum. For k>k, the curve 
labelled e, is €(&, k) defined in (3.129), with @=2ekf. The broken curve 
is &(&, k) with @=2e, where m0. These curves are taken from 
Brueckner and Gammel (1958); they suffer from the error mentioned at 
the beginning of § IV of that paper. 


| Equation (4.116) was evaluated in the singlets-state using the Brueckner— 
Gammel potential given in §1.3. With k;=1-34f-! the inequality was 
satisfied for all choices of e,; with k,= 1-40 f— it was satisfied for the choice 
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(ii) only, and with k,=1-48f— it was not satisfied for either set of é,. As 
pointed out by Mills et al. it is possible that better trial wavefunctions would 
satisfy the inequality in at least some of the other cases. They conclude 
from this work, therefore, that nuclear matter is superfluid. Their 
results show, incidentally, that the effective mass approximation grossly 
overestimates the value of (4.116) compared to the other choices of e,,, and 
is therefore less likely to give a superfluid state. 

Since, as we have seen, the existence of a singularity in the Brueckner K 
matrix already implies superfluidity on the BCS theory, the calculation of 
Emery (1959), referred to in §4.2, supports the conclusion that nuclear 
matter is superfluid. In fact it follows from Emery’s work that a superfluid 
state exists, in the effective mass approximation and at a density corres- 
ponding to k;=1-5f—, provided M*>0-56 M. Mills et al. do not in fact 
obtain superfluidity with these parameters, so their trial functions must be 
inadequate for this case. 

The magnitude of the energy gap caused by the singlet s interaction has 
been estimated by Emery and Sessler (1960). They linearize (4.76) by 
assuming that, in the region around k=; in which p, is not negligible in 
comparison with ¢,, the variation of y, can be ignored, so that it can be 
replaced by Hue = €0 in the demoninator.. From (4.96), 2¢) is the required 
energy gap. Introducing y by (4.110), we can write for the s-wave part of 
(4.76) using operator notation. 


ly)=—@vly) Pegs Sere he (Al D0) 
where 1(? dk|k)<k| 
G=— pea Ne 1S 4.121 
wo Tee eh Boake 
and <r |k) =sin (kr). MAS ee te, Pry wet a soy (4: 122) 
Following Emery and Sessler we define 
|6)=|k,) Gv |4) it gh ree WER SAAICA EE) 
with 1 pdk [| ky <k|— [> dg a: 
Sa eS NE NE Tb Se ge Sat ae 124 
ea) rep eg ane 
Then Ple lx d= Fy lolx <b [eGov | x) 


from (4.123), and 
<b lelxd= — P| eGe |x) 
from (4.120). The difference of the last two equations yields 
ames 7- a), area ihe arob’ 3 
<plvlke> Jo Lex? + eo]! 
Emery and Sessler approximate the integral here by using the effective 
mass approximation, which should not be too bad since the integrand 
is strongly peaked about k=k,, and thus obtain 
8h? k? a ~ he k,* 
ea ape oxP bee one Meee (41 26) 
In evaluating ¢ from (4.123) they ignore the ey in (4.124), which is pre- 
sumably permissible since the integrand is not peaked around k=k,, and 


P.M.S. Z 
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use the form (4.119) for e,, with Cy equal to the binding energy per 
particle and 


They plot the energy gap as a function of M/*/M, using the Brueckner— 
Gammel potential, and with the roughly empirical values k,=1-4f", 
é,.=15-5Mev. Itis a rapidly increasing function, being about 0-1 Mey for 
M*=WM and less than 0-5 x 10-> Mev for M*=0-75 M. At a lower density, 
k,=1-0f-1, with e,, taken as 10mev, they found 2€)=3-5 Mev and 1 Mev 
respectively for the same values of J/*. 

An alternative method of calculating the energy gap is proposed in the 
work of Brueckner, Soda, Anderson and Morel (1960) and Anderson and 
Morel (1960). To establish this method we write (4.76) in the form 


pe SULLY pay os) tebe 46 a ee ea 
where 
k)<k | 
L = L924 bee > . . . . . 4.1: 
> 515 8+|un PP See 
and divide L into two parts 
L=1L,+L, . . . . . . . . (4.129) 


where L, involves states with | ¢, |<de and L, involves states with | «, | >6¢ 
for some (6 e) to be chosen later. Then, from (4.127), 


Lap = Lv (Ly + Lou) 
which gives ; 
Lop = (1—Lyv)* Lav Ly p. 


On using this to eliminate Z,u from (4.127) we obtain 
B=(v+v(1—Lyv) Lev) Lye 
SW GIN Si Pr iN ea Otay ee RUA SSE 


where K=v+v(1—L,v)v 
=V+vUL,0+vLgvL,u+.... Spee ae ie ee 


Anderson and Morel propose to choose de to satisfy (i)5e >, so that p, 
can be ignored in L, and (ii) de sufficiently small for ¢k’ | K|k) to depend 
only on the angle between k and k’ for |, |, |€,-|<8e. If (ii) is satisfied 
then, considering for example just the s-state interaction and putting 
(k’ | K |k) = —v9/Q, (4.130) is equivalent to the simple problem discussed 
in§ 4.5 and, for vy > 0 we have (4.79) 


25e 


2e = —__.,, 
sinh (1/7/09) 


tk aR TO) 
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For v) <0 there is no superfluid solution. ‘Since €y is independent of 6 e, 
(4.79) requires that v, varies rapidly with e for small e and in fact tends to 
infinity asd}e+>0. The solution of the equation for K does not have this 
behaviour however. For example, with a separable potential, it is easy to 
see (cf. the discussion in § 4.2) that «k | K [k) tends to zero as k—>k,. This 
illustrates the breakdown here of (4.79) for small 8e. Probably it is then 
permissible neither to take the matrix elements of K constant nor to 
ignore w in determining L,. 

Brueckner, Soda, Anderson and Morel (1960) make the further approxi- 
mation of equating K with the Brueckner K matrix. Apart from the 
restriction to | ¢,|>8e in K, which is to some extent taken account of in K 
by the method of treating the self-consistency, which introduces a gap in 
the spectrum (see § 4.2), and by the mesh used in the numerical treatment of 
the integral equation, these quantities differ by the inclusion of terms with 
€,<0in K. These correspond roughly to hole-hole scattering terms and, 
although these do not contribute greatly to the total binding energy, they 
may well be significant in individual K -matrix elements. Thus we do not 
regard this approximation as reliable. Brueckner and Gammel found 
{k, | K | k;) to be positive in the singlet s-state and negative in 
the triplet s and singlet d-states. Thus, according to this method, the 
singlet s-state does not lead to superfluidity, a conclusion which is in 
contradiction with that reached above on the basis of more reasonable 
assumptions. For the triplet s and singlet d-states, Brueckner et al. 
calculate the gaps to be 0-0066e and 10-°de respectively. One has little 
idea of the value of 6 e which is appropriate to the particular approximation 
used to calculate K, so it is diffiuclt to give estimates of the energy gap from 
these figures. Brueckner ef al. take de to be the fermi energy (which they 
take as 57-5 Mev) but this is probably too large for the approximations of 
constant matrix elements and K = K to be justified. 

In conclusion then the existence and magnitude of an energy gap in 
nuclear matter is seen to be a delicate question. Moreover the gaps which 
have been found are very small, being at the most about 0-1Mev. The 
mean spacing of single particle levels in a cube of volume equal to that of a 
nucleus with A nucleons is about (30/A) Mev, so that only for the heaviest 
nuclei would such a gap be significant. Emery and Sessler (1960) suggest 
that a larger gap might arise in actual nuclei (see Bohr e¢ al. 1958 for a 
discussion of the empirical facts) because of the lower mean density due to 
the large surface region. Such ideas cannot, however, be taken very 
seriously because the ‘correlation length’ (Bardeen ef al. 1957) is so large. 
In fact the range of momenta involved in the BCS phenomena is of the 
order of M,/k, so that the correlation length, which is the inverse of this, is 


fe Pod 2 eG 
f ~—f 


2Megk, 


where «, is in Mev. For all plausible values of ¢9 this is larger than any 
nucleus. Itisnotsurprising therefore that the pairing effect in nuclei should 
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depend greatly on detailed properties of the particular nucleus considered 
(Beliaev 1958). 

In the volume energy of nuclear matter the BCS contribution per 
particle ismuch smaller than the gap 9, since only asmall number of particles 
near to the fermi surface are involved in the latter. It is unlikely that 
extrapolation from actual nuclei, with subtraction of surface effects, shell 
effects, etc., will ever fix the binding energy of nuclear matter with such a 
precision that such a contribution becomes significant. 


§ 5. CONCLUSION 


The quantitive success of the theory cannot be accurately judged at 
present. On the one hand, the empirical properties of nuclear matter are 
ill-defined (§1.1). This applies even to the binding energy per particle, 
which varies by several Mev between different versions of the mass formula 
(Green 1958, Cameron 1957). On the other hand, although this is not clear 
from the literaturey, the Brueckner-Gammel calculation is sensitive to 
details of the internucleon force not well determined by the scattering 
experiments. In particular the experimental uncertainty in the triplet 
central odd state interaction permits several Mev variation in the calculated 
energy per particle. Thus the numerical results could not, for example, 
rule out the possibility of substantial contributions from many body forcest. 
In forming a conclusion we can only consider whether there exists a theory 
which, were the two body force accurately known, would allow a reliable 
evaluation of the consequences. 

The Jastrow approach is inadequate here. It does not run into the 
singularity troubles of Brueckner theory, and we think it would give 
sensible results if worked out with realistic forces. But these results would 
not be accurate, for the cluster correlations appear to be large. There are 
also other corrections, because of the approximate nature of the basic 
Jastrow wavefunction, which have not been estimated. 

The singularity in Brueckner theory is easily lost in practice, either by 
approximations in the working or, as with Moszkowski and Scott (§ 3.5), 
by not iterating the attractive tail of the potential. The lowest order 
corrections then appear to be small, of order 1 Mev in the binding energy per 
particle at normal density. However, the methods of estimating these 
(§3.6) are very crude. They give, for instance, a much smaller self-con- 
sistency correction then that finally obtained by Brueckner and Gammel. 
These estimates should be repeated more realistically, including in particu- 
lar the tensor force. 

Although the quantitative effect should be very small, the Brueckner 
method will be formally satisfactory only when the singularity is eliminated 
rather than lost. The correct description of the fermi surface is probably 


} See for example the discussion following the report by Brueckner to the 
1960 Kingston conference (Brueckner 1960). 
{ See, in this connection, Derrick et al. (1961). 
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that of Bardeen, Cooper, and Schrieffer, although their theory is 
convincingly established only for a grossly non-local (‘reduced ’) hamilton- 
ian. This account of the long range correlations should be combined with 
the Brueckner account of the short range. This may come about by, for 
example, a coherent formulation of the separate treatment of low and high 
energy virtual excitations (Brueckner et al. 1960), or by a practical version 
of the Green function scheme. 
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